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COAXIAL PLASMA ACCELERATOR STUDY 


ABSTRACT 


This study gives an up-to-date review of the coaxial plasma 
accelerator - introducing regimes, models of regimes, and types of 
accelerators. The theoretical investigation is devoted to the 
axisymmetric-steady-shock-tube model and related numerical procedures 
are developed and verified. The experimental effort is in the 
accelerator with a teflon dielectric. Circuitry for an original version 
of a unipolar pulsed plasma accelerator has been developed as a part of 


this study. 
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CHAPTER I INTRODUCTION 


A coaxial plasma accelerator (CPA in the following text) illustrated 
in a very simplified version by Fig. 1.1 is one member of a large family 
of electromagnetic plasma accelerators. In those devices the acceleration 
of a body of ionized gas occurs by the interaction of currents driven 
through the gas with magnetic fields established either by those currents 
or by external means. For a more complete picture we may refer to Table 1 
taken from Jahn and von Jaskowsky (1968). The subject investigated here 
is marked in the table, namely, we are going to deal with pulsed CPA 
having: static filling of helium, no external magnetic field, internal 
ionization (with small preionization), capacitor-bank type current source, 
direct coupling of discharge to circuit, pure gas and vaporized solid as 
the working fluid, coaxial channel geometry, and with external switch (see 
Fig. 1.1). From Table 1 it is also clear that many conclusions related to 
a CPA are correct also for a parallel rail plasma accelerator, and pinch 
devices, and vice versa. However, generalizations of this nature may be 
misleading. 

Apart from the fact that accelerator phenomena are intrinsically 
interesting, there are some very practical reasons for studying them. First, 
in propulsion applications, this electromagnetic mechanism holds promise 
for providing a combination of high exhaust velocities with high mass flows, 
in comparison with pure electrothermal and electrostatic mechanisms. For 


example, a pulsed thruster is used for station-keeping in the LES-6 
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Fig.1.1 Coaxial plasma accelerator. 
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3 
satellite (Vondra et aZ.1970). These thrusters use an arc discharge across 


a teflon surface;the ablated molecules are then ionized and exhausted at 
3000 m/s. 

The second application area of CPA is in plasma injectors, that is 
production of fast and dense plasmoids. For such purposes the purity of 
the gas to be ionized and accelerated is important. This requires 
investigation of: gas-conditions in the experimental chamber, residual gas, 
fast gas valve operation, gas desorption from the walls and decomposition 
of heavy particles of oil vapors penetrating from the vacuum pumps, etc. 
(Gryzinski, 1968). 

The third application area of CPA or rail plasma accelerators is in 
production of electromagnetically driven shock waves. Here, we require 
the arc to be in the form of a thin impermeable current sheet and to act as 
a magnetic piston preceded by a shock wave. The associated relaxation 
phenomena behind such shocks are studied from a chemical kinetic point of 
view (Chang , 1970). In this way we can obtain much faster shocks than 
with the classical diaphragm version of shock tubes. The incomplete 
knowledge of physical processes involved in these experiments necessarily 
curtails evaluations, especially in this highly promising third area. 

Another application of CPA is the acceleration of small projectiles 
(% 100um) to meteor velocities (% 20 km/s) by exposing them to a stream of 
high density, fast moving plasma (Utah R.D. Co., 1970). A suitable 
insulator between the CPA electrodes could be the source of projectiles. 

These experiments are of current interest, yielding better results 
than "classical" electrostatic accelerators of heavy particles. 

A better understanding of the nature of CPA would be invaluable in 
the designing of such devices. An appreciable amount of experimental data 


and a number of theoretical models are available, but a basis on which to 
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design a particular device having the required operational parameters is 
not available, and each device requires thorough operational testing 
(Sadowski et aZ.1968). 

In comparison with pure electrothermal and electrostatic mechanisms, 
electromagnetic interactions are both technologically and theoretically 
more complex. The situation becomes even more difficult because of the 
transient erosion of the insulator and electrodes. The problem of the 
electron-emission mechanism in an "ordinary" steady arc has still not been 
satisfactorily solved (Hoyaux, 1968). Here we are dealing with transient 
arcs (~ microseconds), moving reasonably fast (~ cm/us) and carrying heavy 
currents (~ 100 kA). Such arcs under certain conditions create a narrow 
Current sheet. Few reliable experimental data exist describing the time 
between application of voltage and beginning of sheet motion. 

This study furnishes an improved model of steady gas flow between the 
current sheet and preceding shock wave (here called the shocked layer model). 
Also, the numerical solution of this model is given for an axisymmetric case 
(ideal gas only), and for a quasi-one-dimensional case (helium with two 
degrees of ionization). It should be pointed out, however, that any model 
describing one particular mode of acceleration is highly limited in use. 

Experimental efforts have been concentrated on two shapes of the teflon 
insulator and two related modes of acceleration have been examined. 

A few comments about the chapter contents are in order. Chapter 2 
is devoted to a description and discussion of physical processes of a real 
but simplified CPA. Chapter 3 discusses the more mathematical treatment 
of shocked layer models and the simplifications necessary for numerical 


solution of those. Some of the difficulties of related numerical 
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calculations are explained in Chapter 4. Calculations in two spatial 
variables are compared with calculations in the literature, and the quasi- 
one-dimensional solution is compared with available experimental data. 
Details of the apparatus (accelerator and diagnostics) are described in 
Chapter 5, which includes an original version of a unipolar pulse gen- 
erator (100 kA,35 psec) with circuit analysis and experimental verification. 
Two completely different modes of acceleration for two teflon insulator 
shapes are described in Chapter 6. Finally, Chapter 7 gives conclusions 


and the author's perspective on the problem. 
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CHAPTER II COAXIAL PLASMA ACCELERATOR 


An incomplete state-of-the-art of CPA is given in this chapter. 
Literature on this subject is extensive and written from various points of 
view although at present, it is impossible to write a unifying theory of 
a CPA. On the other hand, the amount of accumulated information is 
sufficient to give a brief outline of certain typical kinds of behaviour 
(here called regimes) and certain basic types of CPA. 

A review of CPA is arranged in the following manner: 

TYPICAL REGIMES (2.1). A regime is a certain "simple" kind of 
behaviour occuring during a small part of the complete process. 

MODELS OF REGIMES (2.2). A model is usually a set of mathematical 
relations and conditions, and related sets of solutions. These models are 
based on further simplification of physical processes. 

BASIC TYPES OF CPA (2.3). Using the terminology developed in the 
section dealing with regimes, we can give a few basic types of CPA. This 
section is written more generally for pulsed plasma accelerators. 

OTHER PHYSICAL PROCESSES (2.4). This section covers some important 
processes not discussed before. 

The last Section 2.5 gives conclusions regarding our possible work in 


this study. 


2.1 TYPICAL CPA REGIMES 


Each particular accelerator should undergo two or more regimes from 
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the following list. Unfortunately, a real accelerator will perform some 
mixture of these modes and the experimentalist must carefully distinguish 
the importance of each mode involved. This again means that in describing 


these modes, we are dealing with a simplified, idealized accelerator. 


ed | Discharge-initiation Regime 


During this stage the total discharge current should grow from zero 
to, say, one tenth of the maximum value. First, the working medium is 
prepared and ionized. There are three basic arrangements for the discharge 
initiation: (a) A small amount of gas is introduced by a fast-acting valve 
into the evacuated barrel, and ionized. (b) A surface discharge is 
initiated on a nonrefractory insulator, and the ablation product becomes 
@ suitable medium. (c) A metallic foil is clamped between the inner and 
the outer electrode and the interelectrode Space is evacuated. The main 
current pulse first melts and evaporates this foil to produce a suitable 
working medium. 

In all cases a localized discharge is required. Case A is simply 
realized by means of holes for the gas inlet. Case B is improved by a 
suitable ignition electrode. In some cases such triggering is performed 
by a "small" plasma gun. 

The first phase has received very little attention from experimentalists. 
It is an accepted fact that this stage leaves the accelerating gap full of 
unknown preionization, precursors, etc. It is essential to understand 
that this results in serious uncertainty of some parameters in all the 
following regimes. 


Occasionally a bell-shaped distribution of current density several 
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millimeters thick is generated. This process depends on the kind of 
working medium used. Another very important parameter is the current 
rise-time, given by the external circuit design. 

There are no quantitative theories and only rarely are experiments 


devoted specifically to investigation of the resulting current spreading. 


2.1.2 Shock-tube Regime 


Here, we assume that an impermeable thin current sheet has already 
been generated. The situation is much easier if we also assume that the 
current layer is of infinite conductivity. Then the magnetic field (B) 
is non-zero only behind the sheet and its effect can be accounted for by 
the magnetic pressure term, (B°/2u). In other words we have a moving 
magnetic piston. What kind of disturbance will be created is now 
dependent on the piston velocity and the gas medium in front of the 
piston. Usually a dissociating and/or ionizing shock wave is created 


which is separated from the driving layer. 


2.1.3  Plasmoid-acceleration Regime 


For this regime one usually assumes that there is no gas in front of 
or behind the plasmoid. The plasmoid itself is created by the rather broad 
current layer moving along the CPA electrodes. In some models the motion 
is described for the plasmoid center of gravity and the thicker current 
layer is assumed to be localized at the same place, i.e., one dimensional 


models are considered. 
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2.1.4 Fast-plasmoid-creation Regime 


It is known that an important characteristic of the accelerator 
Operation is in the gas-pressure distribution along the system, which 
governs the slow (Section 2.1.3) and fast regimes of accelerator 
operation. The slow regime is characterized by a current sheet of velocity 
10/ m/s and by gas combustion. In the fast regime, most of the plasma 
acceleration (10° m/s) occurs in the density gradient where intense 
polarisation oscillations are produced (Kalmykov et aZ. 1971 and 
Khizhnyak, 1970). 

This topic will not be described further here although much work 


remains to be done. 


2.1.5 Plasma-ejection Regime 


This phase is a transition between the movement inside the channel 
wherein the plasma is created and accelerated, and the ordered streaming 
motion of the gas (plasmoid) outside the accelerator barrel. This regime 
is theoretically very difficult. Experimental effort has been primarily 
concentrated on one special version of this regime, the dense plasma focus 


which is described in the following section. 


2.1.6 Dense-plasma-focus Regime 


In some cases, during the plasma ejection from the end of a CPA, part 
of the stored magnetic energy in the accelerator and external circuit is 


rapidly converted to plasma energy. This process occurs during the sheet's 


yoieraisga, wiz. ta siteiialgam 
Arie ,mazevy2 any PAs mig 


Oe 1S1O5GA TO 


vi PO! BV 70 a es Pah. TY & V0 @Oat'ts 


enétia ahs Po Svea ane ner é 


hannaco. os 
parmsatd2 bosib xh a , oles 
root etat iSty6d tots wa! wage \ silt 
viiveetyy need 


268 24071s ‘fehatntss 


2070? nezg!q Serat) snl , saihgey se a4 


Sore | i elihe a ee ' isu 
< 





} 5 y n a? fe erties As sent AMO er Py i v-~ - 


fab owereyg-esp 500 At 2) ROPERS 
uh (ETS noticag) wolfe 2a) eee 


, Diy HOD é ar Pes Me | A itt . fro sao 


bes! eye Me 


EvaeMOD gap NG 


2 


ft 240300. le\m OL) Tee 
th ila! p22 io ttae beafeg 
(OTR) .Avetiesal 
uorsdsived don, (iy obese sia 


~ | »Sneb sd 09 ener 
SHPDBAr No 77 S8Lar MRT Be 

Wisd Nom 2ngad fet szemyeiies 
os. bte fos bayo ci etiga ig ott Shee 


obizjieot( Stomebl 92) 200. a3. To eRe 
Cisy =] lento 2g 


‘o igh prew-Tetasq2: ene ne bodetnganed 
‘nolizer ootwohlot edd sat; badiewest 2}.hattie 





10 


brief collapse toward the axis. 


10 on 


Typical parameters of the phenomenon are: plasma densities » 10 "cm ~, 
temperatures of a few keV, life-time of 100 nsec, spatial dimensions of a 


few mm? An excellent review is given by Loveberg and Griem (1971). 


2.1.7 Inertial-plasmoid-motion Regime 


Here we assume that a localized amount of plasma (plasmoid) is moving, 
but no external force from any external energy source acts on the plasma. 
During the ejection phase, many "long" extended current loops have been 
established and now these loops are broken. Internal vortex loops may 
still be expected in the plasmoid. 

This regime is usually obtained outside the accelerator barrel and is 
considerably less complex than those regimes dealing with the space 
inside. For the above reason, this phase has received most experimental 
and theoretical effort, primarily because of the plasma-injector applications. 

Let us note that if the current pulse is short and the accelerator 
barrel is long, then such a regime could still be attained inside the 


barrel. 
2.2 MODELS OF SOME CPA REGIMES 


The regimes introduced in Section 2.1 cover a certain important time 
during the acceleration development. We should realize that the real 
physical processes are sometimes more complex. Furthermore, the transition 


from one “simple” regime to another is difficult to study both from the 
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theoretical and experimental point of view. On the other hand, these 


separate regimes are relatively easy to describe using simple models. 


2.2.1 Discharge-initiation Model 


The discharge-initiation regime has been described in Section 2.1.1 
with the conclusion being that there are no suitable quantitative theories 
available. For the Staticcdasepressure type with the perfectly refractory 
insulator, Jahn and von Jaskowsky (1968) give a simple, skin-effect-like 
explanation for a one-dimensional discharge. One can envisage the 
breakdown process beginning with a short-lived uniform glow discharge 
over the entire channel. As the current density rapidly rises, driving 
the discharge through an abnormal glow towards an arc, the associated 
distribution in the time derivative of the magnetic field induces a new 
component of interior electric field, opposing the applied potential and 
preferentially discouraging current flow across the portions of channel 
farthest removed from the external circuit. The current thus intensifies 


most rapidly near the channel entrance. 
2.2.2 Shock-tube Models 


Referring to Section 2.1.2, we know that a current layer acts as a 
magnetic piston, and that a shock wave is driven by this piston. Some 
models of this piston-shock situation are available, and others may be 
created by the extension of some special shock wave models (current sheet 
is then placed on the contact surface position), or by extension of 


current layer models. A short list of such models with brief explanations 
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will follow. GAS-KINETIC MODEL covers such situations where individual 
electron, ion and neutral trajectories may be examined under the pre- 
vailing circumstances of Hall parameters, gyro radii and mean free paths. 
It is known that this approach is very difficult even for only a normal 
plane ionizing shock without current layer and magnetic field (Zel'dovich 
and Raizer, 1966). Collision processes and relaxation-layer discussions 
also belong in this group. | 

MULTI-FLUID MODEL in comparison to the previous case is macroscopic 
as to chemical equilibrium in the region between the driving layer and 
the shock wave. Thermodynamic functions are used to include changes due 
to dissociation, ionization, and excitation of the gas. The presence of a 
magnetic field in this calculated region, depends on the assumed current 
Sheet conductivity, or an external field may be assumed. 

SINGLE-FLUID MODEL with given finite conductivity neglects the multi 
fluid behaviour and chooses a set of equations to represent a single-fluid 
acted on by some source of force. 

GAS-DYNAMIC MODEL treats a current sheet as an impermeable piston 
driving a gas dynamic shock ahead of it (at a somewhat higher velocity) 
into the undisturbed gas. The shock wave is followed by a uniform region 
of hot gas which is bounded by the current sheet. Ordinary Rankine- 
Hugoniot relations for constant specific heat ratio y (in Section 3.2 
called simply shock relations) can be used. The real gas effect is 
expressed by using some effective value of y close to 1.0. 

SNOW-PLOW MODEL is very common in the literature and we may define 
such a model simply as a special case of the previous gas-dynamic model 


with the effective specific ratio being y = 1.0. From the detached shock 
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wave theory we know that this assumption leads to a zero shock-piston 
distance. It is a well knwon fact that the mass of the accelerated layer 
is accumulated from some "initial-breakdown" value. This model will 
appear again in Section 2.2.3 as a limiting case. 

SLUG MODEL is also a special case of the aforementioned gas-dynamic 
model. Here, the gas density in front of the layer is assumed to be zero, 
i.e., no shock wave is created, and no mass accumulation occurs. This 
means that the layer (plasmoid) of constant mass is accelerated along the 
electrodes. This model does not really belong under the heading of shock- 
tube models, but it is the last simplification of these models. 

The previous models could be treated as stationary cases or as time- 
dependent models, in one, quasi-one, or two spatial dimensions. By the 
quasi-one-dimensional model we mean that instead of dependent variables 
Y(r,z), say, we will use variables Y(r,z(r)) = Y"(r). In other words, we 
will get a set of ordinary differential equations instead of a set of 
partial differential equations. 

Let us give a few bibliographical! notes for those models. An excellent 
review of shock-wave phenomena, time dependent one dimensional shock-piston 
models, relaxation effects behind a shock wave including radiation, is 
given by Zel'dovich and Raizer (1966). The problems of establishing 
local thermal equilibrium in a helium plasma created inside the CPA at 
pressures of .5 to 1 Torr, and related relaxation effects are discussed by 
Lie et al.(1967). They conclude that the collisional and the radiative 
processes which determine the population densities of He III, and the 
various bound states of He II and He I are of such a nature that the 


local thermal equilibrium assumption is in fact a satisfactory approx- 
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imation for the conditions of their experiment. This work has been 
extended in a later paper (Lie et az. 1970) in which the temperature was 
measured spectroscopically as well as by laser scattering, where the 
assumption of local thermal equilibrium is not required, and the results 
are in agreement with their previous work. 

The multi-fluid approach with equilibrium chemistry in hydrogen in 
an electromagnetic field is discussed by Taussig (1966). In his 
calculations, no assumption is made concerning that of the contact surface 
being driven by the piston. His results show that solutions including 
chemical reactions are similar to the ideal gas (normal) ionizing shock 
wave solutions. 

The single-fluid model of MHD shock tubes is given by Hoffman (1967). 
His concern was with a diffusing current sheet driving a shock wave, and 
he found that under many conditions, separation is more likely to be 
achieved with a decrease, rather than an increase, in current sheet speed. 
The gas-dynamic model in one dimension and in time is given by Morse (1962). 
The energy source consists of one charged capacitor and a constant specific 
heat ratio is assumed. Approximate closed-form solutions, for "small" 
and "long" time, are also given. 

A quasi-one-dimensional, stationary, gas-dynamic model of the CPA has 
been developed by Pert (1968). He also gives some analytical solutions to 
this problem. The same author (Pert, 1969) extended his previous model 
for real gas effects in hydrogen when equilibrium single ionization is 
included in the energy equation. Butler et az. (1969) give a model 
similar to Pert (1969) but use onlya threshold approximation for a "general" 


Single ionization. 
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Two-dimensional time dependent gas dynamic (ideal gas) calculations 
are also described by Butler et al. (1969). The snow-plow model is very 
simple and still found to be adequate for some shock-tube investigations 
(Chang and Kofoed, 1968). The "old" slug model remains in use too. 
Vondra et aZ. (1970) discussed the teflon pulsed thruster by using this 
model. It will be shown in this Study that in such cases the current 
is distributed along the electrodes, and that the existence of a Slug of 


a plasma is unlikely. 
2.2.3 Plasmoid-acceleration Model 


Following is a description of the time dependent linearized model of 
the CPA. Such an accelerator consists, for example, of parallel-plate 
electrodes with inductance L, per unit length, and a capacitor Cy 
charged to voltage Ue (see Fig. 2.1). A very narrow current layer of 
mass m, is created at position x = 0 and time t = 0. The layer (plasmoid) 
is ideally conducting and is accelerated by the electrodynamic force Za 


and decelerated by the drag force F The insulator is assumed to be 


drag’ 
perfectly refractory, i.e., no ablation products and hence no currents 
near this insulator for t > 0. 

The momentum conservation law and two circuit equations are usually 


written in the following form (Artsimovich, 1963): 
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dts pgrey-y 


dt =e) (2.3) 


plas 
where v = dx/dt is the plasmoid velocity, L = Ly + Lix is total inductance 


of the circuit and V represents some plasmoid voltage drop. 


plas 


The model description starts from the electrodynamic force Fad 


created by the magnetic pressure Be /2u (u is permeability), acting on the 


left-hand-side plasmoid surface of area Dd ai eo 
F_, = bd B¢/2u (2.4) 
ed ° : 


Assuming that b >> d, using (ITT Co., 1969), one can write the following 


equations for the magnetic field, B, which is inside the accelerator: 

B = ul/2 (2303 
and the following for the inductance per unit length, 

L, = udb, (2.6) 
and finally from Eqs. (2.4), (2.5) and (2.6) the electrodynamic force is 


L, 1° (2a) 


Some basic assumptions regarding the mass rate processes must be 


included: 
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Fig.2.1 Sketch of the linearized CPA, 
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Fig.2.2 Types of solutions to plasmoid-acceleration model. 
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dt = SUCKS Vel oe.) (2.8) 


as well as those concerned with the nature of the drag force F For 


drag’ 
the simplest version, the SLUG MODEL, one obtains the following equations 


dm _ 4 
Cie aera Oe (2.9) 


Some mass accumulation is contained in the SNOW-PLOW MODEL, with the 


following assumption: 


m=m, + bdox, or dm/dt = bdov (2.10) 


Farag ot 

As previously mentioned, the hypersonic current layer drives a shock 
wave, with the gas in the layer-shock region (shocked layer) being very 
hot. Hence, the gas is ionized and dissociated, with an effective specific 
heat ratio y = 1.0. Some authors (for example, Kolesnikov and Stolovich, 
1970; K-S henceforth) neglect this type of mass accumulation. 

The undisturbed gas effect is partially covered by the drag-force 
term 


2 


drag 


This is simply a dynamic pressure VS times the plasmoid surface bd. On 


the other hand, combining Eqs. (2.1) and (2.10) one gets 


oo 
{ 
7 


a 4 (8.$) (on Ld ¢ e 
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abd d@x a Za 
(m, 2 Sy ge ppiaal eta (2.12) 
i.e., one has derived an effective drag term and with the addition due to 
K-S, the mass accumulation is included. 


More sophisticated mechanisms of the mass rate processes and drag 


forces have been developed during the 1962-1972 period, and are summarized 


by K-S: 
dm _ 2 2 
qe = 7aym - apne + a,|1| crag aes ies) 
F = bv + bovm + ball |v (2.14) 
drag ] (a 3 ° : 


The aym term accounts for particle diffusion, charge-exchange and electron 
attachment; the agin term describes mass decrease due to particle re- 
combination; the a |1| and Swit terms respectively describe erosion mass 
addition due to ion bombardment, and electrode melting due to Joule 

heating. The drag force terms denote the Poneaince biv represents friction 
between the moving plasma and electrodes and bovm + b.|1| expresses friction 
due to mass transfer. The coefficients a5 b. are purely empirical, others 
could be estimated from simplified mechanisms (K-S). 

The usual manner of solving Eqs. (2.1) to (2.3) is to transform those 
equations into some dimensionless form, to choose initial conditions and 
typical parameters, and then to solve that set of ordinary differential 
equations. Figure 2.2 exhibits such solutions. For details see K-S and 


many other authors. 
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An important conclusion from these figures is that the resulting 
solution [m(t), v(t), x(t), I(t)] is an extremely complicated function of 
the parameters (a, 4-++2b3.0.V, 0,5...) with the theoretical estimations 
at best being very rough. Parameter estimation from experimental measure- 
ment is difficult and is more dependent on the assumptions than on the true 
physical picture. For example, Sarah (1969) using the snow-plow model 
assumes only electrode erosion. Similarly, different experimental 
conditions dictate a vaniety of explanations for limited plasmoid 


velocities. 
2.2.4 Current-sheet Models 


Tne present topic is concerned with a different approach to plasmoid 
acceleration, which has had an equally long history. These more 
sophisticated models make experimental verification more difficult. The 
problem is sketched in Fig. 2.3. 

Before proceeding further let us introduce some classifications of 
the idealized situation. The regions #1 and #2 of Fig. 2.3, are represented 
by corresponding gas density Pa» gas velocity Vi» gas pressure p., magnetic 
field B., electric field E., and electric conductivity O.. The following 
classification of Burton (1968) shown in Table 2 is not generally accepted, 
but on the other hand, such a classification clarifies other papers in this 
field. 

Models of this type are available in MHD single and two-fluid 
approximations, for collision dominated electron motion, and for ionization 


through the current sheet. A critical review of those models is given by 
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Pert (1970 a,b). 


Leaving out mathematical detatls of the above models, let us give 
the resulting description of the sheet acceleration with the collision 
dominated electron motion (Khiznyak and Kalmykov, 1967). One assumes 
that the current sheet is already created. The electrodynamic force 
jxB acts directly on the electron component and accelerates the current 
front along the electrodes. Neutral gas of region #1 is ionized by the 
current front and increases the plasmoid mass. The ions, dragged by the 
established polarization front attain one half the velocity of the 
polarization front, with most of the plasma being dragged behind the current 
front. The polarization front creates longitudinal electrostatic waves in 
the plasma concentrating the plasma in potential wells. Since the axial 
electric field has an oscillatory character, the current forms a series 
of closed loops within the plasma. The total discharge current I is equal 
only to the current flowing in the front and tail of the plasmoid. Electro- 
dynamic acceleration (Section 2.2.3) occurs only when the energy of these 
electrostatic oscillations is small in comparison to the energy of the 
directed plasmoid motion. These electrostatic oscillations are intensified 
by driving the current sheet through the gas density gradient (Khizhnyak, 
1970). The latter case results in creation of a fast moving plasmoid. 

This model does not forecast any mass accumulation in front of the 
fast current layer. Conclusions of other models (Hoffman, 1967) are 
similar. In recent papers such separations have been experimentally 
observed, and a missing mass-accumulation mechanism has been determined: 
resonant charge-exchange collisions of neutrals with ions (Pert, 1970c 
and Val'kov et al. 1970). 


A rather simplified, but very common description of the current sheet 
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Table 2 Classification of steady MHD discontinuities. 
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Fig.2.3 Diagram of current sheet motion. 
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Fig.2.4 Sketch of insulator ablation model for one time instant. 
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assumes that ionization occurs at the front of the sheet. Then the 
electron Hall parameter A ge: 1 and the ion Hall parameter wits < ] 


’ 


whereby, all current is carried by electrons. Thus 
J, 7 en, , (2.15) 


1.e., the electron current density We is due to the ExB drift of the 
electrons in the crossed polarization field E and magnetic field B (Jahn 


and von Jaskowsky, 1968; Pert, 1970c; and Valkov et aZ. 1970). 
2.3 BASIC PULSED-PLASMA-ACCELERATOR TYPES 


As has been noted in the introduction, some basic features are common 
for several pulsed plasma accelerators (PPA in the following text). This 
means that after introducing the typical regimes and related models, we 
can more easily describe some basic types of PPA. 

It is useful to roughly divide these types into two groups. The first 
covers those used to supply a gas or plasma with energy or momentum in 
a relatively short time, and the inertial-plasmoid motion regime is then 
established. This group includes the various forms of T-tubes, and of 
conical, theta, and Z pinch devices. The second group is comprised of 
those systems in which the accelerated plasma is continuously supplied 
with energy and momentum so that a steady state can be achieved with the 
right kind of energy storage (current pulse generator). Typical members 
of the second group are the CPA and the parallel accelerator, sometimes 


called coaxial and rail gun respectively. The first group is reviewed by 
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Muntenbruch (1969). This study is concentrated on the second group. 

One could mention that the division of the PPA in these two groups is 

not the best, because the same short CPA, for example, belongs to the second 
group for the plasmoid acceleration regime, but it belongs to the first 


group for the inertial-plasmoid-motion regime. 


2.3.1 Pulsed-gas Type 


Here, the space between relatively long (coaxial) electrodes is 
evacuated. Then a puff of gas is introduced by a fast-acting valve into 
the barrel space, through the holes in the central electrode near the 
insulator. The gas has approximately 100usec to spread out.Now the 
pulse of current is introduced, i.e., the discharge initiation regime is 
performed and a current layer is created in the neighbourhood of the 
admitting holes. The remaining gas forms an ionized gas cloud and the 
current layer drives a shock wave through this cloud, i.e., the shock-tube 
regime is obtained. When the current layer, or better the current front, 
reaches the end of the cloud the plasmoid acceleration and ejection regimes 
should follow. 

The most common example of this type is the plasma injector, designed 
for generation of hot plasma for thermonuclear research. Different 
modifications of the CPA are used for this purpose. One advanced form has 
the outer and inner electrodes made from rods, i.e., these electrodes are 
"transparent" for the particles rotating in the magnetic field. The 
technical difficulties connected with the construction of such a device 


are explained by Gryzinski et az. (1968). 
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@.3.2 Dielectric-erosion Type 


There are many known geometrical versions of this type. The feature 
in common is a nonrefractory dielectric placed between two metallic 
electrodes. An auxiliary ignition electrode starts a surface breakdown 
followed by a rapidly increased insulator erosion (ablation). This is 
caused by surface arc discharges, supplied from the external current pulse 
generator. This is known as a dielectric-erosion kind of the discharge- 
initiation regime of a PPA. 

The next step depends on the current pulse shape, on the energy 
supply generator, and on the shape of the metallic electrodes. For example, 
10 cm long rail plate electrodes are used for pulsed solid thrusters 


(Vondra et az. 1970). 


2.3.3  Electrode-erosion Type 


The PPA of this kind is also called the end type PPA, because plasma 
is formed from the eroded material from the ends of the two metallic 
electrodes, which are separated by a heat-resistant dielectric. The 


erosion mechanisms are of current interest (Khizhnyak et aZ. 1971). 


2.3.4 Shock-tube Type 


Basically, for this we assume that the accelerator has two long 
electrodes, the space Beeieeh is filled by a certain gas of known pressure 
and the current layer (magnetic piston) is created during some suitable 
discharge-intitiation regime. The required current pulse is long and 


unipolar, if possible. Now, the conditions should be such that a well 
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separated shock wave is driven at a constant speed by the current layer. 
In the literature, the electromagnetic (electric, magnetic) shock tube 
is a device where the inertially moving plasmoid drives a shock wave 


outside of the accelerator. 


2.3.5  Static-gas-pressure Type 


We can easily show that this type is similar to one or more of the 
previous types depending on electrode geometry, insulator type, gas 
pressure, and current puise shape. This name has been introduced because 


this type is very common in experimental research of the PPA. 


2.4 OTHER PHYSICAL PROCESSES 


In the previous part of Chapter 2 we introduced some regimes in CPA 
behaviour and escribed a few current models for those regimes. Later, we 
introduced basic types of PPA. Further details of the physical processes 
involved follow below. 

PRECURSOR MECHANISMS. In a realistic model of PPA the effect of 
precursors should be taken into account. From this point of view, the 
experimental results and theoretical conclusions for classical pressure 
driven shock waves, given by Weyman (1969), are of interest. For one set 
of conditions he found that electron diffusion is the major factor in the 
production of precursors. After improving the vacuum conditions, photo- 
jonization was the dominant factor in precursor production. 


Lubin (1966) studied the precursor waves in an electromagnetic T-tube 
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(hydrogen and argon, 1 Torr); his conclusions were in close agreement with 
the nonlinear ionizing electromagnetic wave model. The same author regards 
the following mechanisms as unimportant: electron diffusion from the main 
shock wave or from the discharge region, soft X-rays from electron bombard- 
ment of the anode, the ionizing shock waves. Many authors report experiment- 
al results with photo excitation and ionization as the most possible 
explanation (Fearn, 1968). 

RELAXATION EFFECTS. In many models of electromagnetic accelerators, 
the relaxation effects are neglected. Lie et al. (1970) calculated 
ionization relaxation times in a homogeneous, transient plasma. Since 
collisional ionization occurs mainly via intermediate excited states, the 
jonization relaxation time is determined by the slowest process in the 
excitation chain. For single ionization of helium initially at a pressure 


8 to 1077 sec. For double 


3 


of 1 Torr, they observed a relaxation time of 10 


17 


ionization they observed 3 x 10°” sec for n > 10 “can ~ to several 


e 
microseconds for Sis eds 10/6 cm7>, Comparing these times with plasma 
lifetime (107®sec) the second case is quite high. Relaxation behind 
strong shock waves in a monatomic gas and in multi-atomic gases is 
discussed by Drawin (1971). Ionization by precursors and by intial 
breakdown of the "undisturbed" gas in the accelerator is highly important, 
but usually unknown. When unspecified impurity atoms are present, the 
jonization process may change. Argyropoules et az. (1969) made a study of 
the coupled distribution of current density, electron temperature, and 
plasma composition in non-equilibrium plasma flows, and how these depend 
on relaxation effects arising from finite reaction rates and electron- 
energy convection. 


DISCHARGE PINCHING. The spatially periodic structure of a heavy 
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surface coating, by adsorbed gas or different oxides, will be stabilized 
by those preliminary discharges. The fact that a boundary, which exists 
between the plasma and electrodes, weakens present theories is generally 
accepted (Hoyaux, 1968). 

Conditioning of the glass insulator plate for the vacuum spark gap 
is described by Mather in Lovberg and Griem (1971) as the formation of 
numerous metallic sites none of which overlap on the glass surface. These 
metallic sites disperse uniformly over the surface and form a microscopically 
graded insulator, which can sustain large electric fields without break- 
down. Conditioning mechanism in the CPA may be similar. 

PLASMA FLOW BEHIND THE CURRENT SHEET. Toroidal paravortices are 
formed in the wake of a moving current sheet in a coaxial plasma gun. The 
axial Hall eiectric field, proportional to the axial driving force, falls 
off as re from the central electrode. The vector sum of this Hall field, 
and a radial electric field which falls off as rol, gives rise to a net 
electric field of over-increasing tilt as a function of radius. The 
electric particle drift caused by this field pattern produces the 
necessary swirl in the wake. That is, considerations of conservation of 
mass and momentum flow in the CPA make it apparent that plasma vortices 
are produced behind the current sheet. The above conclusions of vortices 
measurements are due to Bostick et az. (1965). A similar current vortex 
structure, behind the current sheet in a parallel plate rail gun, is 
reported by Pert (1970c). This vortex must be distinguished from the gas 
flow "bubbles", which are periodically formed and detached from the 
discharge in the rail and coaxial guns. For the rail gun, this latter 


phenomenon is limited to a narrow pressure region (Pert, 1970d). Drift 
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velocity of the plasma, behind the current sheet and the vortices is 
given by the relation 
+> 
U, = 3 (2.16) 

where E is the radial electric field and B is the azimuthal magnetic 
field. This velocity, for a short part of the discharge current duration, 
may be either in opposite direction to the current sheet velocity 
(Zolototrubov et aZ. 1965), or in the same direction and of a greater 
magnitude than the current sheet velocity (thus called a fast plasma 
component). The motion of a single charged particle inside a CPA with 
static electric and magnetic fields, E and B, is considered more 
generally by Kulinski (1966). Some of his conclusions may be used for 
certain cases of flow behind the current sheet. 

INSULATOR ABLATION. Theoretical investigation of this effect on the 
performance of magnetic piston shock tubes is given by Workman (1965). 
A related model is sketched in Fig. 2.4. A simple correlation-parameter 
model of Workman involves the magnetic field, undisturbed gas density, and 
ionization energy of the insulator vapor. The above model permits 
determination of the loss in a given experiment. 

Lie (1970), reporting results related to refractory and non-refractory 
insulators, claims Workman's theory is insufficient because the electrode- 
radius-ratio is not taken into account. Still no improved version of the 


insulator-ablation model is available. 


2.5 CONCLUSIONS 


Regimes and types of the pulsed CPA have been given in Chapter 2. 
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3] 
Study of the CPA is more than ten years Old, yet this quick and in- 
complete review shows us that the experimental and theoretical invest- 
igation of the CPA is far from being satisfactory. In the experimental 
field, the dielectric erosion type of accelerator is of increased 


interest, so two shapes of teflon insulator are investigated in this 


study. 
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CHAPTER III THEORETICAL DESCRIPTION OF THE SHOCK-TUBE REGIME 


Previous investigators have studied the shock tube in various 
approximations of physical processes involved, from gas-kinetic to 
snow-plow models. Different numbers of independent variables have been 
used, from the simple time dependent model in one Space variable, Y(z(t)), 
to the time dependent model in two space variables, Yir2et) eric 
object of this chapter is to analyze the Steady (stationary) gas flow in 
the region between the idealized curved driving current sheet and the 
shock wave. 

The problem, formulated in Section 3.1, gives the conservation laws 
in integral form. Such a form can be attacked via different levels of 
numerical approximation. The form also allows both gas dynamic and 
chemical-equilibrium expressions for some dependent variables. 

Equilibrium conditions, for one special case behind the plane shock 
wave, are given in Section 3.2. Two axisymmetric versions of the resulting 
model are then developed. The two-dimensional, and the quasi-one- 


dimensional models are described in Section 3.4. 
3.1] FORMULATION OF THE STEADY SHOCK-TUBE MODEL 
Consider the plasma acceleration between two coaxial electrodes with 


inner radius ce and outer radius >> fo? To carry out the analysis, we 


further make the following assumptions, not strictly rigorously formulated. 
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The situation is illustrated in Fig. 3.1. 
AS1. The current sheet is thin, ideally conducting, already created, 
Stable and axially symmetrical, and moving at a constant speed a in 
the -z axis direction. 
AS2. The current sheet acts as an impermeable (magnetic) piston driving 
a (detached) shock wave into the undisturbed gas. The gas in front of 
this shock wave is not preheated (not ionized) and the gas behind the 
shock wave is in some specified equilibrium. The relaxation layer between 
those two states can be neglected. 
AS3. The conductivity of the gas in the region between the current sheet 
and the shock front (shocked gas, or shocked layer) is high enough to en- 
sure that the driving current I is confined within this thin sheet. 
AS4. The layers adjacent to the electrodes are thin in comparison with 
the distance between the electrodes and the effect of these layers can be 
ignored. 
AS5. An external generator provides a current pulse such that the current 
I confined within the current sheet is constant during the time of 
investigation. 


AS6. The magnetic field, B, behind the current sheet is 
B(r) = pl/2nr (33) 


where py is the permeability, I the driving current, and r the radius 
(independent variable). Because the current sheet is ideally conducting 
and thin, the effect of this magnetic field can be reduced to some magnetic 


pressure Pmag? acting on the sheet: 
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OUTER ELECTRODE p.h,p 
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1 
8 Cae A SS Pmag mee 
UNIONIZED Qn} CONTROL VOLUME 


. GAS r 
r 
Vo 
Po c 


SHOCKED LAYER 
SHOCK WAVE ——> 


IMPERMEABLE CURRENT INNER ELECTRODE 
SHEET 


AXIS OF SYMMETRY 


Fig.3.1 Sketch of part of a long CPA with the control volume just between the current sheet 
and the shock wave;aotation used in the quasi-one-dimensional model in the cylindrical 
coordinates (r,z,). 


(r+ 5 r)(b+.5 b) 





Fig.3.2 The control volume of steady shocked layer model.The surface of this volume is divided 
into six parts with unit normal vectors fy,...,n6 3(r,2,¢) are cylindrical coordinates,v, is the 
axial velocity of the undisturbed gas relative to the current sheet,A and @ are defined by Fig.3.1 
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3D 
Sephe Bo/2u = w12/8ner? (3.2) 
AS7. The density of the shocked gas is sufficiently high that the 
following integral form of the conservation of mass, momentum and energy 


can be used for any control volume inside a shocked layer (Becker, 1968): 


fovpds = 0, 


s 


Jie, nods = 00, 


S 


foth +h ds =o. (3.3) 

2 n 

S 

The examples of the control volume are shown in Figs. 3.1 and 3.2. 
In Eq. (3.3), S is the control volume surface, V is the velocity, fi is 
the outward normal unit vector, NE ss V.U is the normal component of the 
velocity, p is the mass density, p is the pressure, and h is the 
enthalpy per unit mass. 


AS8. Two state equations are provided, for example, in the form: 


p= plo,ie 
h = h(g.T) . (3.4) 


where T is the equilibrium temperature of the gas. 
AS9. The problem is axisymmetrical and using Fig. 3.3, the boundary 


conditions are chosen: 
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Fig.3.3 Sketch of the boundaries. 
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Fig.3.4 Calculated number densities n,n,,and ny behind a plane shock wave of velocity Vo in 
helium of 1 Torr pressure.The critical velocities v,, and v., were introduced by Butler et al. 
(1969).*,values of n calculated by Horton and Menard(1969 ; all excitation levels included. 


+,values of n calculated by Fucks and Artmann(1963);wrong algorithm. 





290N 77342 RITUO 








FOMNTIZET AsV 


a ind 


37 

aan eat) 
#2: oblique shock relation (local approximation) 
#3: arbitrary (transparent outer electrode) (oss 

#4: ecRU se bas Prag 

The oblique shock relations are discussed in Section 3.2. 
AS1O. The effect of viscosity and thermal conductivity is neglected. 
The form of Eq. (3.3) can be rewritten in an approximate form. With 
reference to Fig. 3.2, the surface S of the control volume consists of 


six elements 6S;. Now, the approximate form is 


ee p.V,.6S, = 0, 
a _ 
ele WV5PG%yq * NyP4)8S; = 9, (3.6) 
ne C 
Re eis 


where the subscript i means average value over the surface element 6S. 


3.2 EQUILIBRIUM CONDITIONS BEHIND A STEADY SHOCK WAVE 


Boundary conditions for the shock wave can be obtained by choosing 
a suitable control volume (see Fig. 3.2), placing the boundary between 
the surfaces 2 and 4, and then limiting the thickness A to zero. 
Assuming that the total control volume surface S = }éS; is small then 


the resulting relations are common oblique shock relations (Becker, 1968): 
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Vot Vot ; 
PoYon ~ Pavan » 

(Be7) 
Povon * Po ~ Pavan + Po » 


where Not is the total (stagnation) enthalpy and Vi is the maximum 
velocity for expansion into vacuum. The last two values are used for 


normalization in some dimensionless models. 
3.2.1 Calorically Ideal Gas 

Assuming that the gas is calorically ideal then the specific heat 
ratio is usually 1.67, 1.40 and 1.33 for monatomic, diatomic,and poly- 
atomic gas respectively. The velocity of sound, a = (yp/o)'/% defines 
the Mach number 

M=v/a or My = igh (3.8) 

Finally, the equations of state are 


Sey sal Ree z (3.9) 


where R = 0.8313 x 10° msec”! °K! is the universal gas constant and m 
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39 
is the dimensionless molecular weight. Now, the boundary conditions, 


Eq. (3.7), can be rewritten in this form 
Po/P6 =a F4(yoM, 56) ry 
Po/Po SS Foy M8) ’ (3.10) 
Sue f3(y M26) . 

Details are readily available in Becker (1968). 


3.2.2 Monatomic Gas 


Equations (3.4) have a special form for equilibrium ionization of a 
monatomic gas. An excellent review of the equilibrium and nonequilibrium 
plasmas, including an equilibrium-condition discussion is due to Drawin 
(1971). Only those relations used for further calculations are presented 
in the following work. 

Assuming the down-stream gas satisfies local thermal equilibrium 
conditions, then for a monatomic gas the so-called mass action law has a 


familiar Saha form: 


int s/c 
Mt] Me r+] 2P 4, (T) (2nm,kT) ER rel - AE 4] 
a Ts Sy (Tn) 5s iNtyowde, ao ee exp (- T ) 
r oe (T) h 


(3a) 


where subscripts or superscripts r and (r+1) are related to ionization 


stages r and (r+1), pint is partition function, ER rt] is the ideal 
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is the lowering of ionization energy, Pp is the Debye radius, and other 


ionization energy of the particle in ionization stage r,AF 


symbols have their usual meanings. 

Reference to Fig. 3.4, where our calculations for helium are presented, 
makes clear that the differences between the simplified approach 
ia) =*1* pe: = ae oe = 1) and the more sophisticated (all excitation 
levels included) is not of basic importance in the calculation. 


Introducing the number density of heavy particles: n = No Felaete ls cee 


] 2 
and the average mass of heavy particles m, the mass density p, pressure Ds 


and specific enthalpy h are now given by: 


CHES, 
where r = (e,0,1,2,..) is the summation index. In the enthalpy relation, 
the corrections due to excitated states are usually neglected (Fucks and 
Artmann, 1963; F-A hereafter). 

Now, having the required equations for pressure and enthalpy, Eqs. 
(3.12), we can combine then with Eqs. (3.7) and (3.11), and the problem 
of determination of the local thermal equilibrium behind a (locally) oblique 
shock wave can be solved. 

The solution of the above equations has been perfected for the case of 
a normal plane shock wave in helium at pressures of 0.1 to 10 Torr. Figures 


3.4 to 3.5 show results of the numerical solution for the case close to the 
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42 
experimental situation studied here. The numerical constants are those 
of F-A, but the results are different (see Fig. 3.4). The reason for 
the discrepancies is found in the numerical solution of F-A (Artmann, 
1963). Other authors have also discussed difficulties of this kind 
(Taussig, 1966; Lapworth, 1969). The main numerical problem is placed 
in the transition region between the first and second ionization. 
Solutions of this problem are available for many gases (Horton e¢ al. 
1969) but only for a limited shock velocity region; and in helium only 
for single ionization. 

The second purpose of the calculation illustrated in Figs. 3.4 to 
3.6 was to have a simple and verified subroutine for the equilibrium 
condition calculation. Such a subroutine can be used for the two-and 


quasi-one-dimensional model solutions. 
3.3 TWO DIMENSIONAL MODEL 


The problem formulated in Section 3.1 is known as a blunt body 
problem with given distribution of the pressure on the body (in this case 
the current sheet). Equations (3.5) and (3.4) are one version of the two- 
dimensional model. The form of Eqs. (3.4) depend on the assumptions 
regarding the gas. The form of Eqs. (3.6), however, is closely related 
to the method of solution. The integral form is best only for some 
numerical algorithms (Richtmyer, 1968). Usually, suitable curvilinear 
coordinates are introduced and Eqs. (3.6) are transformed into a set of 
(partial) differential equations. The transformation is done By choosing 


a coordinate oriented control volume in the limit of zero volume. 
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Fig.3.7 Diagram of the axisymmetric flow system showing the variables and notation used.In 


(B),the variables (r,8@) are replaced by variables (y,s).The region of interest is represented 
by points Proeee Py ° 
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Fig.4.1 Comparison of u(s,y) calculated by LS-method with those calculated by Belotserkovskii 
et al. (1966) ;sphere,ideal gas,M> = 10 ,y=L.4. 





Fig.4.2 Comparison of -v(s,y) calculated by LS-method with those calculated by Belotserkovskii 
et al. (1966) ;sphere,ideal gas,M, = OR View la 4 i 
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Fig.4.3 Comparison of p(s,y) calculated by LS-method with those 
: al 
et al. (1966) ;sphere,ideal gas,M.= 10 ,y= 1.4. se calculated by Belotserkovskii 


For spherical coordinates, the results of this transformation are 


found to be (Bird et al. 1960): 
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po(v += + — 


where v and u are now the components of the velocity in the r and 6 
directions, see Fig. 3.7. These equations do not contain the transformation 
of the energy equation. The reason for this is that, using an appropriate 
assumption, one can replace the differential equation by a known integral, 


for example, by the Bernoulli equation (Pert, 1969): 


2) = const (sai4) 


Details of the numerical solution are discussed in Section 4.1. 
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3.4 QUASI-ONE-DIMENSIONAL MODEL 


The previous two-dimensional model will be shown to require a lengthy 
numerical solution. To simplify the problem, we will assume that all 
dependent variables of the shocked layer are functions of one space 
variable, r, say. In other words, the shocked layer is assumed to be 


thin. Referring to Fig. 3.1 where the cylindrical coordinates (r,z,0) 
are used, this also means that the control volumes of Fig.) 3-2 nas the 


surface element 8S, on the current sheet and 5S5 On the shock wave. 

One can imagine a simple algorithm, solving numerically Eqs. (3.4) 
and (3.6) for a series of points placed in the middle of the control 
volumes between the sheet and shock. More common, however, are methods 
for solution of ordinary differential equations. Using a procedure similar 
to that in Section 3.3, one can write Eqs. (3.6) for the control volume 
of Figs. 3.1 and 3.2, and in the limit of 6rs0, the following set of 


Ordinary differential equations is obtained: 





fed 
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Po (hy + 7 Volo Sin [o(h + x V hare ohlpes (3.18) 


where the subscript o refers to undisturbed gas, p = nm is gas density, 


v is gas velocity, A is shocked layer thickness, h is specific enthalpy, 
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(r,6) are coordinates defined by Fig. 3.13; p 1s pressure, p is magnetic 


mag 
pressure, and m is the mass of heavy particles. 

Using the assumptions that AO 3) = 0 (stagnation line near electrode) 
and that hy << we, the mass and energy conservation laws can be integrated 


to give 


(3719) 


Note that the last equation is again the Bernoulli equation (3.14) in 
changed notation. 
The z coordinate of the middle of the shocked layer can be related 


to the (r,6@) coordinates (see Fig. 3.1) by the relation 


= tane (3720) 
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Now, Our model contains Eqs. (3.16), (3.17), (3.19), and (3.20). The 
equilibrium conditions of the gas should be added, i.e., for the ideal- 
gas by Eqs. (3.9), or for equilibrium ionization of a monatomic gas by 
Bas.. (3211). and (3.12). 

Unfortunately, the model is still incomplete (more unknowns than 
relations) due to the quasi-one-dimensional assumption. This can be 
overcome, however, by making an additional assumption concerning the 
solution. Accordingly, the closing relation is taken to be that used by 


Butler et aZ. (1969): 
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p = 5 (0,5 cos“e + Prag? (3e71) 


Here, the magnetic pressure Pad is given by Eq. (3.1). In other 
words, the pressure in the layer is approximated by an average of dynamic 
pressure on the front and magnetic pressure on the back. 

Numerical solution of Eqs. (3.11), (3.12), (3.16), (3.17) and (3.19) 


to (3.21) is discussed in Section 4.2. 
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CHAPTER IV NUMERICAL TECHNIQUES 


The steady shock-tube models, described in the previous chapter, 
involve some specific problems investigated in this chapter. General 
formulation of the model transformed our problem into a numerically 
difficult group of problems, called the blunt body problem (Hayes and 
Probstein, 1968). 

Calculation of supersonic compressible fluid flow around blunt 
bodies has received considerable attention in recent years. Several 
successful numerical calculations have been made using techniques like: 
(i) Finite difference method (Gentry et al. 1966; Amsden, 1966) 

(ii) Method of integral relations (Belotserkovskii, 1966). 
(iii) Pade's expansion method (Van Tuyl, 1970) 

A common feature of those methods is that they are very complicated 
and sensitive to variation of some parameters. For example, to make the 
method of finite differences stable, one has to use certain "artificial 
diffusion", introduced only for such a reason (Van Leer, 1969). 

In Section 4.1 is described a variational method, the discrete non- 
linear least squares (NLS) method, applied to this problem. Such an 
approach has not appeared in the literature, to the author's knowledge. 

The NLS method in our modification is written with two independent 
variables and involves a two level minimization procedure. Two level 
minimization places strict requirements on the convergence of the method. 


In the case discussed here it is done by a two level non-linear least 
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Squares procedure. Non-linear approaches are based on a "sensible" 

first guess to the solution, followed by iterative improvements of the 
weighting coefficients or functions (in the case of more than one 
condition/equation) to obtain the "exact" solution. A general description 
of the variational approach may be found in Mikhlin (1971). 

The advantages of this method are: 

a) In the case of trans-sonic flow the method is stable through the sonic 
line; the discontinuity at this point is a problem in many Eulerian 
or Lagrangian techniques, but may also be overcome by using a mixed 
coordinate system such as the P.I.C code (Amsden, 1966). 

b) The solution cf the problem for the entire flow region is described 
by a small number of parameters (10 to 50), so that the method has 
only small memory requirements. | 

c) Any known physical properties of the solution may be easily included 
in the form of the approximating functions, e.g. if the functions are 
even or odd. | 

d) The problem can be easily inverted (e.g. given shock wave, pressure 
distribution etc., calculate shape of body etc). 

In Section 4.2, the same problem is simplified by using a quasi-one- 
dimensional assumption, i.e., the set of ordinary differential equations 
is to be solved subject to the initial conditions. The Runge-Kutta method 
(R-K) is used in combination with the solution of a simultaneous set of 
nonlinear equations (equilibrium conditions). 

Initial conditions have to be found by using a suitable assumption 
and then replacing the differential equation by an auxiliary set of additional 
nonlinear equations. Solution of the extended set gives the required initial 
conditions. The last procedure physically means that the flow has a 


stagnation line near the central electrode and, to find the gas parameters 
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and the shock-sheet distance, we need to make an assumption regarding the 
resulting flow in the region near the stagnation line. 

The NLS solution is done only for very simplified cases and compared 
with known solutions. The R-K solution is compared with existing 


experimental results in the literature. 


4.1 LEAST-SQUARES SOLUTION TO THE TWO-DIMENSIONAL MODEL 


In this section, one numerical experiment is described attempting to 
solve the steady shocked layer model in two spational variables. It is 
assumed: (i) the current sheet is spherical and the central electrode 
radius is zero, (ii) gas has constant specific heat ratio y. 

A variational method, the NLS method is applied to this problem. The 
objective of this particular numerical experiment is to provide "adjustable" 
simplification to fluid-dynamic programs of some laboratories (say, of Los 
Alamos; Butler, 1969). 

The basic idea of a variational approach is the transformation of a 
boundary value problem into a problem of finding the minimum of certain 
functionals. The previous assumptions transformed our problem into 
hypersonic flow around a sphere, whose solution is known (Belotserkovskii, 
1966). 

The model consists of Eqs. (3.9), (3.13), and (3.14), for example. 
These equations can be rewritten in the form of three differential equations 


for unknowns Y(u,v,p) with independent variables x = (r,6): 
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fy(xsY) = (ov). + (ou), (2v + u cotan @) 
f,(x,Y) = rp,./e TAIN not URES u2 = 0 
Ie ep, Pat UU ete ryur tony = 0 (Aoi) 


where partial derivatives are denoted by subscripts and where 


Pemorey hamock Neg teh (4.2) 


The quantities u, Vv, p, o and h are two components of velocity, pressure, 


sup ve iz 5 


mass density and enthalpy respectively, normalized to v max?fo” max 


max? i 


Py and nis respectively. The subscript o is related to free stream 
values, velocity Veet related to total enthalpy Not: Functions fi: fos 
f, have been introduced for further use in the variational reformulation 
of the problem. 

The sketch of the flow system and notation used is in Fig. 3.7. Polar 


coordinates (r,e) are replaced by the "body-shaped" coordinates (y,s) 


according to the relations: 


y = (r = reqpy)/ (shock ~ “Bony? 


8 
b. 2 2,1/2 
S = fcr rpopy * (draopy/de) aude (4.3) 
fe) 


where s is the length of the rpopy (°) curve. Assuming that reopy = i, 
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52 
in our case, the latter equation is reduced to s = 6. The purpose of 
this transformation is that the approximate functions in the next section 
are simpler in (y,s) variables. 

The boundary conditions, Eqs. (3.5), can be rewritten into the 


following form: 


v= 0 for r = rpopy: U = 0 for 6 = 0 (4.4) 
and 
Youock'Po? Nos Vos o) for r = reugey - or 
(4.5) 


SHOCK ~ 


where Youock denotes the value just behind the shock wave, the angle o 
is shown in Fig. 3.7. Note that it is not possible in this particular 
case to write the boundary conditions for 6 = Oras which are required if 
the finite difference approach is to be employed. 

A great deal of theoretical work has been done on the problem of 
finding sufficiently wide classes of functions in order that some functionals 
actually attain their minimum (Mikhlin, 1971). This section presents two 
basic modifications of that "theoretical" approach: 
i) The class of problems considered will be the one for which the 

uM 


av] 
approximate form of the solution, Y(x,A) = (Yi Yoo ¥ ),; is known. 


For example, let 
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avers 2 n 2 4 2 
Pw ¥3 = Mig t Agg® + Apy@” + (Aga + AygO” + Ang® Jy + (Ang + Anes 
+ Aog8 )¥* 
uta of eount ised 
shock » 4] * Age” + AS 
(4.6) 


ii) In this analysis of the detached shock wave problem, the position 
and shape of one part of the boundary (shock wave) is given implicitly by 
the Hugoniot relation, Eqs. (4.5). 

We have already mentioned the method of integral relations, which 
has three versions (Scheme I, II, III). Our proposed method differs from 
the third version in the form of the functions given by Eqs. (4.6) and.in 


the method by which the coefficients A and A” are determined. 
4.1.1 Variational Formulation 


Let A” and the points x be given. The solution of the boundary value 
problem, Eqs. (4.1) to (4.5), is given by the set A for which the functional 


given by 
av) 
F(A) = ye we eee Y) (4.7) 


is minimal. The determination of weighting factors w, will be discussed 
later. 


A second functional 


F (A ) = F(A(A )) (4.8) 
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is introduced for obtaining the solution for A. under the assumption that 
the correct shock shape and position is Obtained when F" is minimized. 
The usual method (Mikhlin, 1971), for minimizing F and F involves 


the solution of equations 


oF 
FARh i = 0, Q = 25 »N 
3A, 
* 
LP ee ‘< 
ang OFeike ty2 5 (4.9) 


Where the derivatives are replaced by numerical differentiation. If N 

is set to 27 in Eqs. (4.6), the boundary conditions, Eqs. (4.4) and (4.5) 
may be used to reduce N to 15. Providing a reasonable first guess for A 
and A is available, the solution to Eqs. (4.7) and (4.8) may be obatined 
by the two level minimization of F and Ea using some generalized Newton 
method (e.g., Newton-Kantorovich method; Mikhlin, 1971). The convergence 
of the above method, however, was found to be slow. An alternate technique 
was used to accomplish the minimization, and such a technique is described 
in Section 4.1.2. 

The weighting factors We in Eqs. (4.7) are determined by using the 
assumption that the first guesses for A and A are good for at least small 
Ae (i.e., near axis) as follows: 

i) For the grid in Fig. 3.7 values for (xs 5¥) are evaluated with Ae = Ae/10. 


ay) Ww. are then determined from the solution of the equations 


-" 7) : ° , | ; 
4 Aen betanulevs svn 4 on vent: san nr bing 
7 


® 


{ ; . 
7 a 


Hd Gps S113, 70 nota ioe wy ig ee Picard: 
a hk 
i. i re ; 
= ve 





55 


4.1.2 Minimization of the Function of N-Variables 


Minimization of the functionals F(A) and F(A) of Section 474.1 
is described. Both the following algorithms have been developed for 
this study. 
i) Minimization of F(A). A more rapid convergence can be obtained if a 
better approximation to the minimized functional is available. This may 
be achieved by linearizing the functions FLLX, » V(x, 5A)], ae i Wee 
1 = 1,2,...,M; rather than linearizing the functions of aF/3A, 3 
2 = 1,2,...,N; in Eq. (4.7). This is performed by using the weighted non- 
linear least squares approach to solve the system of 3M equations 


av) 


FLEX; » Y(x; sA)] athe (4.11) 


dhe. ae Ages ++ An. This is accomplished by iterative improvements of a 


] 
given first guess ae through the solution of a linear least squares problem 
in the neighbourhood of Sb 


ii) Minimization of aCee The minimization of the second functional can 


be performed as follows: 
7 zs * * * * 
a) in the neighbourhood of the first guess (Al, Ay, Ax) = (xy Xos 
* * 2 
X3)» N > 10 values of F are evaluated, 
b) Ft is then approximated by the polynomial 


2 2 2 | 
Fi wek(ajes a]X7 + aoXy + agxXq + AgXyXo + ALXIX, + AGXOXg + AX) 


+a + AgX3 + a4 (4.12) 


8X2 


Pee epee 
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given by the soiution of 
av) 
3F/ ax, so OR SS A ere fo (4.13) 


if a, < 0 for some i€(1,...,6), the next guess can be determined using 


Only first derivatives. 


4017.3. -Resikl ts and Comparison 


The method described in the previous parts of Section 4.1 have been 
used to solve cases with specific heat ratio y = 1.4 (% air) and with Mach 


number M = 10. The chosen weighting functions were 
W (y 8) = 0.15, Wo(y26) = 0.61, W3(y,8) = 0.77 


i.e., only constants without iterative improvement, and three iterations of 


F(A) have been used, M = 2, N= 15. The resulting shock wave was 


2 4 


= U.175 +70.) le. + 0,044.5 (4.14) 


"SHOCK 
The numerical results are compared with those obtained by the integral 
relation method, version III the so called method of transcendental 
equations; Belotserkovskii, 1966). Figures 4.1 to 4.3 indicate that the 
agreement is satisfactory. 
The following conclusions can be drawn from the above analysis: 
i) Detached shock wave problems can be solved using the least squares 


approach. 
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Fig.4.4 Initial condition calculation;A(r, + é6r) versus driving current; quasi-one- 
dimensional model. 
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Fig.4.5 Ideal steady-state current sheet velocity versus driving current;points US ae a 
are of the same meaning as those in Figs. 4.7 to 4.14. 
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ii) The calculation time is typically 2 minutes on IBM 360/67 for the 
case of a sphere in anideal gas with the number of representing points 
M = 24 (see Fig. 3.7); number of A coefficients, in Eqs. (4.6), N = 15; 
and with three iterations for three A coefficients. This time may be 
improved by reducing the number of Xa s Ass Ay used, or else using a 
faster minimization procedure. 
iii) The solution is relatively insensitive to the values of weighting 
factors Whee but for 

Wy) = My = Wy = (1/3) 144 
the procedure is unstable. Note that our polynomial approximations, Eas. 
(4.6), are written for the whole region. Using the local polynomial 
approximation (spline functions, Greville, 1969) one may improve the 
stability. 
iv) The resulting algorithm, only roughly described in Section 4.1 and 
by Jurak,Pert,and Capjack (1972) is rather complicated. The flow ina 
CPA, under more realistic conditions, will be attacked by using the quasi- 


one-dimensional simplification in Section 4.2. 
4.2 RUNGE-KUTTA SOLUTION OF THE QUASI-ONE-DIMENSTONAL MODEL 


One can easily see that the solution of the two-dimensional problem 
could be rather lengthy. The quasi-one-dimensional version, hoever, makes 
the situation much easier (Jurak and Offenberger, 1972). The model from 


Section 3.1 can be rewritten in form: 


) SRT YO HOFTE 


) 302 AT. TOY: 
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de _ (a cose)/re - a/(r2cose) 2 99 
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se (a sine)/ré - (a v)/ (vere) - A dp/dr 





Fa es, 9) leprae ep aes eer cress ra et a arnt eam (4.15) 
dr es 
2 2 
are POO 0 
ae tan 6, A a om 


ny(n + 2 9) - 2c(n - ties ng)T/ “exp(-E,/T) = 0 


ae 
No (ny + 2no) - C nT / exp(-E,/T) = 0 (4.16) 
a cos“6 ] 
where the pressure p = DAlgocks =) and the total number density 
4 i 
) 


n=p/m=n +n, + ns. 
o/ 3 1 + 
Comparing with Section 3.2 where mks units were used, here cgs are 


employed with the exception of temperature T and ionization potentials 


E, and E, which are in eV, and current I which is in Ash = iG 107!4 


21 3-3/2 for helium. 


erg/eV and c = 6.04 x 10° cm “eV. 
| Our problem is to solve a simultaneous set of three ordinary 
differential equations dY/dr = Fi(rsYsP), where Y = (6,v,z), and 

P = (To 5m) sn.) is given by the solution of the set of nonlinear equations 


FA(P) = 0. 
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Bee lenin Lialecondi Fions 


The initial value of v is given according to assumption AS9: 


i 
fa) 


vite) 
It is also assumed that 


a(r,) : 


i} 
(=e; 


Zr,) = Alro)/2 ; 
where z is the coordinate of the middle of the shocked layer.Integration 


can not start from point r=r, , because 1/v(ro) = » in Eqs.(4.15). 


The following approximations to the first two differential equations were 


made 


deere bye! Or 
dr ér 6r 
REND Ve 
dr ~*~ é6r ér 


which is "correct" only near r = beets: Choosing sufficiently small or 
(say 6r = 0.01), the complete set of model equations is now a set of 
nonlinear equations. Figure 4.4 shows one curve of this solution giving 


shocked layer thickness A(r.) as a function of the driving current I. 


4.2.2 Integration of the Model 


Knowing the initial conditions, we can use common Runge-Kutta procedures 
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for integration of the set of first order differential equations. The 
only difficulty was that for each integration substep the set of nonlinear 
equations had to be solved. An algorithm for nonlinear equations has 
been developed as a modification of the Seidel method (Ortega e¢ al. 1970). 
Note that the problem of instability of the common Newton-Raphson method 
for equilibrium calculations (Taussig, 1966) has been overcome by using 
different (equivalent) forms of the same set of nonlinear equaitons for 
different regions of ionization. This algorithm was first verified by 
solving the problem of equilibrium ionization behind the plane shock wave 
in helaum (Figs. 3.4-to 3.6). 

Numerical stability of the latter procedure was verified by using a 
"very bad guess"; the incorrect solution of the same problem is also 
illustrated in Fig. 3.4, and the iterative improvement of the guess value 


to the correct value is shown in Fig. 4.6. 


4.2.3 Results and Comparison 


Figures 4.7 to 4.14 show the results of the numerical solution of the 
model for helium at 1 Torr over a range of velocities which covers no ion- 
ization to doubly ionized. One important result of this model is (Fig. 
4.13) that the shocked layer thickness does not go to zero at some radius 
as found by Butler et aZ. (1969). The model is very simple (less than 1 
minute of computer time), so that extension to cases with different 
assumptions concerning ionization equilibrium or energy conservation are 
easily accomplished. 

Our CPA, described in Chapter 5, was not operating in the required 


shock-tube regime, but there are available suitable experimental results 
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Bo 


Ld 


TOTAL NUMBER DENSITY (1017 cm’) 





0 100 200 300 400 


NUMBER OF ITERATIONS 


Fig.4.6 Iterative improvement of the guess value.The resulting value is one point 
on the n = n(v,) curve of Fig.3.4 Se = 4.5 cm/sec). 


cm ?) 


TOTAL DENSITY (10!7 





RADIUS (cm) 


Fig.4.7 Calculated number density n = Nye + Nyt + nytt versus radius ;quasi-one-dimensional 
model.See also Fig.4.5. 
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Fig.4.8 Calculated number density ny = Th versus radius ;quasi-one-dimensional-model. 
See also Fig.4.5 
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Fig.4.9 Calculated number density ng = nygt+ versus radius;quasi-one-dimensional model. 
See also Fig.4.5 
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COORDINATE z (cm) 
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Fig.4.10 Calculated equilibrium temperature T versus radius; 
quasi-one-dimensional model.See also Fig.4.5 


1,2,3,4,5 


RADIUS (cm) 


Fig.4.11 Calculated z coordinate of the middle of the shocked layer versus radius; 
quasi-one-dimensional model. See also Fig.4.5 
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Fig.4.12 Calculated incidence angle @ of the middle of the shocked layer versus radius; 
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quasi-one-dimensional model. See also Fig.4.5 
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Fig.4.13 Calculated shocked layer thickness versus radius;quasi-one-dimensional model. 
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See also Fig.4.5 
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Fig.4.14 Calculated gas velocity along the current layer versus radius; 


quasi-one-dimensional model.See also Fig.4.5 


of Lie et az. (1967). In this work the authors reported spectroscopic 
measurement of electron temperature and density for plasmas produced in 
a CPA filled with helium to a static pressure of .5 to 1 Torr. Luminous 
front velocities of ~ 6 cm/us were attained. Both spectroscopic 
measurements and our calculations depend on the local thermal equilibrium 
assumption. Lie et al. concluded, after careful discussion, that for the 
conditions of the experiment (and calculations) such an assumption is 
in fact a satisfactory approximation. 

Figures 4.15 and 4.16 give comparison of the experimental data with 
the numerical solution. The experimental electron number density, for 


negative inner electrode, shows reasonable agreement with calculated 1 


n, + 2ny. Note that the reported measurements were done with a plastic 
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NEGATIVE 
INNER 
ELECTRODE 


MEASURED CALCULATED 


POSITIVE 
INNER 
ELECTRODE 


ELECTRON DENSITY (1027 cm 3) 





RADIUS/r, 


Fig.4.15 Calculated and experimental electron number densities;v. = 6 cm/Usec ; experimental 
data of Lie et al.(1967). 2 


(Acrylic) insulator and copper electrodes (r = 1.5 cm, r = 2.5 cm), 


inner outer 


and performed just behind the end of the outer electrode. The inner 
electrode was longer to prevent collapse of the sheet. This means that 

the electron density inside could be in even better agreement. The 
Calculated and measured equilibrium temperature, for negative inner 
electrode, are in excellent correlation. It is of interest that the 
comparison is done for velocity Pe = 6 cm/usec where both degrees of 
ionization should be taken into account. There is no model in the available 
literature suitable for such a situation. The model of Pert (1969) 

requires some modifications and a numerical solution while the model of 


Butler et aZ, (1969) gives only a very simplified ionization approximation. 
















jwivensy _ 
: : {hii a : ; 
|;  200nt02A . 
} | é 
, ; : = 
AY i bat - 
OsTAIUIIAD Of \ ‘GIRWRAIN —- 
y “/ = 
/ z j 
ff = ‘ fa 
a f ' r1204 & 
YA _ Finest — 
ys ™hATS3 a 2 : 
L Wt a 
' 
; a _ _ hes =o ee ee 6 
f : v , 
Puri 


lasooetshexe ; 15 \oo ¢ v srerdh te0ee woe ia sceoce, -oghe Le wsaleatad ‘Pisa 
4 f*a¥l).ia Jn-all Seo 


» (ms €.% = edo’. <7? {= seus? 2a soe teyq > bre-rosefuent (>) Irak) 
yenat SAT .aboutaoto vsdwe at? tobe odd bniitad 3 2ut dsirrot 9g: ns, 

dad seyen 2taT .deadé-2n2 Yo Geqal foo demvenc at noynot 256 Sato 

— , ae Pe Sees uriansd news 

; Fat ae ups berwesen bas beds sfusts 

deen Yo af 2% peateveas fas} feowe nt ois (bbe 


ce 
towt rE A 


















68 


CALCULATED 


X 


NEGATIVE 
MEASURED INNER 


ELECTRODE 


TEMPERATURE (eV) 





RADIUS/r,, 


Fig.4.16 Caiculated and experimental equilibrium temperature; Vea 6 cm/sec , experi- 
mental data of Lie et al.(1967). 
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forming network using two trigatrons for unipolar pulse generation. 


CHAPTER V APPARATUS 


69 


The experimental part of this study is devoted to an investigation 


of the acceleration regime of a CPA provided by a teflon insulator. 


block diagram of the apparatus used for this study is in Fig. 5.1. 


TRIGGER 
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CURRENT PULSE GENERATOR 
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Fig.5.1 Schematic layout of the CPA apparatus. 


Section 5.1 deals with the mechanical design of important parts of the 


apparatus. 


Section 5.2 describes the design and operation of a current pulse 


the first trigatron has external triggering. The resulting current 


A 


Table 3 shows some details and parameters. 


Only 


— ee 
; 





MULT SAY 
WhO?) 


D 





CUTARAGSR: Vo ATTIAHD 


ob 2f ybute side Yo deg Lesnmiiroqes ont 


vd bebivewg ATG 5 oe eafasy nortedeleses GAP Te 


wie =i) “ot boro eus cage GH To mempeth Soon 


ure AU AINSO a) AY Fea? *, 


2 






= 





70 


pulse is approximately 100 kA of 35 usec duration and has a shorter 
rise than fall-time. Calculated characteristics show good agreement 
with measured values. 

Section 5.3 gives some details of diagnostics, infrared interfero- 


metry and luminosity, visible luminosity using photodiode, and magnetic 


probes. 


Table 3 Description of the apparatus. 









COPPER,OUTER DIAMETER = 1.58 cm 
INNER ELECTRODE LENGTH = 60 cm oibela syne 
OUTER ELECTRODE COPPER,INNER DIAMETER = 5.0 cm Fin sa 







LENGTH = 60 cm 






INSULATOR | 





TEFLON ;D, = 1.59 cm, D> = 4.99 cm 


S 


INSULATOR 11 TEFLON;a = 6 mm, b = 4.3 mm 








10°? Torr,EDWARDS ROTARY PUMP ED 150 





BASE PRESSURE 


STATIC GAS PRESSURE 






HELIUM AT PRESSURE 0.1 TO 10 Torr 





SPARKS OVER THREE AUXILIARY TEFLON SURFACES; Fig.5.] 
CAPACITORS: <9 3x" 1OcnFo/ 1S kV Fig.6.1 


CURRENT UNIPOLAR;100 kA , 35 usec Frigesal 
PULSE GENERATOR TEN CAPACITORS; 14.5 uF / 20 kV , 50 nH Fig.5.4 
0 TO 10 kV de Fi 
| Gesu 
THYRATRON; 300 V PULSE GENERATOR 
TRIGGERING CIRCUITRY 300 V DELAY UNITS Fig.5.1 
10 kV PULSE AMPLIFIER 


TRIGATRON;10 cm DELRON CHAMBER,AIR 1 atm 
SWITCHING ELEMENTS 2.5 cm COPPER ELECTRODES 


PREIONI ZATION 


2 mm COPPER TRIG. PIN PLACED 
INSIDE ONE ELECTRODE 
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o.1 MECHANICAL DESIGN OF THE ACCELERATOR 


Arrangement of the apparatus is illustrated in Fig. 5.2, and the 
accelerator electrodes are shown in Fig. 5.3. The first version, GUN 1 
had an outer electrode made of brass which was "transparent" to particles 


rotating in a magnetic field and for diagnostics. However, because of 
poor breakdown characteristics a second, GUN 2, was designed. This 


version had windows on one side only and the mechanical Stability was 
improved. Finally, the third version, GUN 3, had two windows, but on 
both sides, i.e., the outer electrode allowed active optical diagnostics 
(external light source). The pyrex-glass cross (Fig. 5.3) with plate- 
glass windows was placed in the same axial position. For infrared 


diagnostics a NaCl window replaced the plate-glass. 


VOLTAGE CURRENT ACCELERATOR INSIDE 
PROBE PROBE OF 7.6 cm PYREX 


GLASS TUBES 
EXTENSION | 


COAXIAL. [[o10( Ae ean ed iy. Se ee : 
LINE 


—— ae ae 
' 









TEN CAPACITORS 
14.5 uF/ 20 kV 


ROTARY 
PUMP 


RESISTOR Ro 
(2.5 cm STAINLESS 
STRIP LINE) 





Fig.5.2 Sketch of the mechanical set-up.Notation according Fig.5.4 








2 , 
ae 
ee Le 
ee, - 
-_- ie Inks ¢ 7 
be O. Hay ¢, Pe: Mm Li ‘ 


| | 











id 
7 


vat ms weR fi 


stents a = 


0 ne bat | ; 


7 
nofeisy in 


ee 
aVO*Oh! 
a 


: bated 5 
ard TA, 


72 


TRANSPARENT" OUTER ELECTRODE 








GUN 1 
OUTER ELECTRODE 
GUNG 2° 7 .. - =. = 
t Qo 4 
We INNER 
ELECTRODE 
INSULATOR 1] SPURIOUS 
BREAKDOWN 
INSULATOR 


/ 
y7 TRIGGER PIN TERMINALS if 





Fig.5.3 Three versions of the CPA. 
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9.2 ELECTRICAL DESIGN OF THE ACCELERATOR 


For these studies a plasma accelerator with long electrodes is 
required, i.€., a high current (100 kA) pulse generator with known rise- 
time, long duration, unipolar behaviour and fall-time greater than rise- 
time, Unipolarity is important for creating only one current sheet per 
shot which avoids interaction between sheets. In addition, observation 
of electrode erosion is considerably simplified (Section 7.1) and one has 
the possibility of performing integrated spectral measurements. 

Such a long unipolar current pulse has been generated by means of a 
charged LC delay line and three switches. This type of device, providing 
currents of a few kiloamperes, was built in CERN (Caris and Williams, 1969) 
for a different application. Each of the switches had a pulse generator 
for triggering and delay unit for time adjustment. For plasma accelerators, 
however ,use of a matched delay line is complicated by the inductive 
character of the load (approximately few hundred nanohenries and few 
milliohms). The characteristic impedance Lo of the delay line is of the 


order of 100 milliohms, which means that a matched delay line results ina 


considerably lower current. 

A simple solution, for currents of a few amperes, can be obtained 
with a charged delay line approximately matched on the other end by a 
resistor in series with a reversed diode. This is not feasible, however, 
for high voltage (10 kV), high current (100 kA) applications as considered 
here. Another more complicated solution uses two sets of capacitors and 
switches (Mikoshiba et al. 1967), or one capacitor and time-varying 
resistance wire (Cheng,1969). 


The method discussed in this section is a mixture and extension of 
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the first and the third schemes mentioned (Jurak and Offenberger, 1972a). 
A charged delay line (ten capacitors) is discharged into the plasma 
accelerator and in the middle of the first half-wave a second switch is 
opened on the opposite side of the same line into approximately one half 
the characteristic impedance. This second switch is triggered by the 


wave generated on the active side of the line. 





Fig.5.4 Circuit diagram of the pulse generator for a unipolar pulsed plasma accelerator. 
The plasma gun is approximated by an inductance Ly sR* = 1M2 ,C* = 10 nF »R, = 10 kf. 


0.2.1 Circuit Analysis 


Figure 5.4 shows a simplified diagram of the device. The circuit 


can be described by the following equations: 
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and two relations: 


A pes 


t 10. (LC) Zemn(L (cyl? (5.2) 


where the moving current sheet in the plasma accelerator is approximated 


by 
Vt (5.3) 


i.e., by an inductance varying linearly with time. 


The first switch provides initial conditions of 


Ie =, 1 ees (5.4) 


and the second switch is simplified by means of the condition that the 
; * 
derivative of Ih, is zero up to a time t (triggering time of the second 


switch). 


5.2.2 Experimental Results 


These equations were numerically integrated using the Runge-Kutta 


method, the results of which are shown in Fig. 5.5, for Ro >> La: The 
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CURRENT 1, (kA) 


0 


VOLTAGE Uyq (kV) 





0 20 40 60 80 100 


TIME (usec) 


Fig.5.5 Calculated I,and Uj for parameters L)v = 5.6mH/sec,L = 196 nH ,L, = 98. nH,L, = 111 nd, 


R= 8 m2,R, >> 2, = 118m2,C = 14.5 UF,V, =10 kV (in this case t* has no influence on the 
results). 


calculated current and voltage waveforms show good agreement with the 
experimentally measured waveforms of Figs. 5.7 and 5.8. 

In addition, the dependence of currents Ty Ios Iy3 On shorting 
impedance Ro is illustrated in Fig. 5.6, which shows that the ideal value 
of Z,/Ro for optimum unipolar pulse generation is ~ 1.8 for our given 
experimental parameters. 

The experimentally observed current waveforms when a nearly correct 
value of Z/ Ro was chosen are shown in Fig. 5.8. Though the discharge 
current waveform is not quite optimum, since the exact value of Z/Ro 
was not used, the result is nearly unipolar, which is adequate for some 


situations. A slightly different choise of parameters would have 
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provided better matching and consequently improved current pulse. 
The results show in any event that the technique is simple and able 
to generate long unipolar high current pulses for plasma accelerators and 


can be easily modified for different applications. 
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CHARACTERISTIC IMPEDANCE / Ro 


Fig.5.6 Variation of currents Iyi »Im2 ,and 143 with Z,/R2,for experimental parameters as above. 
These currents are defined in Fig.5.5. 





Fig.5.7 Oscillogram of the current I and the voltage U,93R,>>Z, Ve = +8 kV,1I,:100 kA/div, 
Uy9:5 kV/div,10 usec/div. 
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Fig.5.8 Oscillogram of the CULrentse Ls sand: 1 
3k, v .8,V, = -8 kV,Iy: vy. 
ois ; ai ay om 1 11> 2 Z ipa (Sh Pay 8 k »1,:50 kA/di tay is not 


9.3 DIAGNOSTICS 


This section describes the details of the diagnostic methods used in 


this study. 


9.3.1 Current and Voltage of the Discharge 


Figure 5.9 shows a schematic of this measurement. The CPA is preceded 
by a 15 cm long coaxial line, where the high voltage probe (Tektronix 
P6015) has a special connector for safe and high frequency voltage pick-up. 
A Rogowski-coil, in classical toroidal form, is also built into the 
coaxial extension line. The signal from this coil is integrated to give 
a value proportional to the total discharge current. 

Voltage is measured with a Tektronix probe. The current probe is 
calibrated by replacing the accelerator with a stainless steel strip-line 
resistor R of unknown value. One capacitor (C = 14.5 uF) from the 
capacitor bank is used and charged to Us = ] kV and discharged into R. 


The voltage U, (see Fig. 5.9) is proportional to the total current I which 


is given by 
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where T = 2L/R, w = fey ales ps A U fol. Determining the value of k 


is straightforward. 
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50 2 VOLTAGE PROBE 
Fig.5.9 Schematic and typical oscillogram of the current and high-voltage probe.Five 
capacitors charged to 10 kV.Current:50 kA/div,voltage:2 kV/div.Horisontal scale:5 usec/div. 


For the current probe calibration the accelerator was replaced by strip-line resistor 
and only one capacitor was used. 


5.3.2 Luminous Front 


The visible light from one spot to the inter-electrode gap is chosen 
by two distant iris diaphragms (3 mm diam) and picked-up by a semiconductor 
photo-diode (Fig. 5.10). Using different windows (holes) in the outer 
electrode, the diode axial position is changed for successive shots, and 
the resulting multiexposure of oscillograms (Fig. 5.10) gives the time 
deyelopment of the luminous front. This measurement yields high shot to 
shot reproducibility. Figure 5.11 illustrates modification of the previous 


method, where the latter assumption (reproducibility) is not required for 
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the velocity determination. 


Using NaCl windows and a Ge:Au detector cooled by liquid nitrogen, 


we can also measure the infrared luminosity. 


ACCELERATOR 


* CURRENT PROBE 





MIRROR Dor SPLITTER 
| 
| 
v 


hpa 4205 


2 kQ 
PHOTODIODE 
+20 V TO SCOPE 
1 uF . 


Fig.5.10 Luminous layer position measurement; experimental set-up and typical oscillograms; 


V, = -10 kV,I:100kA/div, photodiode:10 v/div,horisontal scale: 10 usec/div.One photodiode and 
six exposures were used. 
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Fig.5.11 Luminous layer velocity measuremen;experimental set-up and typical oscillogran, 
I: 100 kA/div,photodiode:0.2 V/div,horisontal scale:5 usec/div. 
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9.3.3 Current Distribution 


The velocity of the current front can be measured by one or more 
non-integrated magnetic probes (10 turns, 1 mm? cross-section) placed in 
4 mm pyrex-glass tubing (Fig. 5.12). The current distribution along the 
electrodes, or the current sheet profile, can be measured in the same 
Way by using two opposite coils, 1 cm apart (Fig. 5.13). The latter 
version is sometimes called a differential magnetic probe and its theory 
is given by Belan et al. (1970). Note that even the classical magnetic 
probe is of interest for studying fast discharges. Pachner (1971) has 


recently completed a doctoral dissertation on this subject. 


4 mm PYREX TUBING 
P 


CURRENT PROBE - 
100 kA/div 


TWO COILS 
Al = 10 cm 


CURRENT PROBE 
50 kA/div 


FOUR COILS 
Al = 1.25 cm 





Fig.5.12 Current sheet velocity measurement;schematic and typical oscillograms for two 
and four coils nonintegrated magnetic probes (10 V/div);horisontal scale:5 usec/div. 
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Fig.5.13 Schematic and typical oscillograms of the differential magnetic probe; I:l00kA/div, 
magnetic probe:2 V/div,horisontal scale:2 usec/div,axial position:6 cm. 


5.3.4 Infrared Diagnostics 


Three types of infrared measurements will be discussed. CO, LASER 
INTERFEROMETRY employs a Mach-Zehnder interferometer where the wavelength, 
4, used is 10.6 um. The i r region requires special material for windows 
and beam splitters (usually NaCl and Ge). Since details are readily 
available for pulsed as well as c.w. measurements (Offenberger and Kerr, 
1972 and Offenberger, Kerr and Smy, 1972), only pertinent relations are 
given. 


The integrated optical path 
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SU a firelsrat (5.6) 


0 


to produce one fringe shift at = 10.6 uM WS 92.1 x 10!6 cm72, A 
necessary condition for this measurement is that w of the radiation be 
greater than the plasma frequency Wy i.e., the upper limit of electron 
density for A = 10.6 um is ccalol? cm7 3, 

The experimental arrangement is shown in Fig. 5.14) Splitters and 
Windows are made of NaCl, and the mirrors are Silvered-glass plates. 

An idealized output from the detector for a linear Change in electron 
density Ne in time is illustrated in Fig. 5.15. 

A schematic diagram of the nitrogen cooled Ge-Au detector is given 
in Fig. 5.16, where the cooled photoconductive detector has a nominal 
resistance of * 100 kn. Time response is limited mainly by stray 
Capacitance since the intrinsic time constant is ~w~ 20 nsec. The detector 
itself has a capacitance of 10 pF. The arrangement in Fig. 5.16 provides 
for three outputs with different source impedances. Here, the time 
response is limited by the connecting (coaxial) cable between output and 
oscilloscope. 

Alignment of the interferometer requires that the narrow beams 
(diam ~ 3 mm) create beat signals on the oscilloscope, i.e., the small 
area of the detector (% 2 mm) has to pick-up only one fringe at a time. 
C0. LASER ABSORPTION is based on 10.6 ym radiation absorption via inverse 
bremsstrahlung (Offenberger and Kerr, 1972). The transmission T is 
related to the absorption coefficient a by 


L 


T = exp(~ f a(s)d2) (5.7) 
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Fig.5.14 Mach-Zehnder interferometer arrangement. 
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Fig.5.15 Evaluation of the interferogram for a linear electron density change and path 
Al ~ 4.2 cm.One fringe shift corresponds to 2An, . 
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Fig.5.16 Schematic diagram of the Ge-Au detector. 
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where £ is the path length, and a (cm7!) is related to the electron 


number density No (cm7) and electron temperature Pavey) for A = 10.6 um 
z eB oe3/2 
& =" 1282x710 pale (5.8) 


It has been assumed that hv << kT, where hv is the photon energy and kT 
is related to the mean thermal energy per particle. Note that for temp- 
eratures of a few eV, measurable absorption occurs for nN, 29 Xx jon eccnit> 
Optical components, C0, laser, and detector are as in the previous method. 
BREMSSTRAHLUNG RADIATION EMISSION is a method opposite to the previous 
one (Managadze et al. 1968). According to theory, the emitted power, 


3 sec!) of ir bremsstrahlung radiation in a frequency 


Qdv (erg cm 
interval dv is related to electron density No (cm7) , ion density 


n (cm7) , and electron temperature ie (°K) by the relation 


where g is the Gaunt factor (=3) and z is the number of elementary charges 
of the ion. Again it is assumed that hv << kT and that the light has only 
a continuum spectrum in the region of interest (narrow ir filters are 
to be used). 

This method can be used as an absolute (with necessary calibration) 
Or relative measurement. For the relative measurements, the detector 


Signal, Uy» resulting from radiation from volume é(vol), transmitted 
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through a filter with bandwidth é6v, and assuming that ne Nae Ss oiven 


a | 
by Eq. (5.9) 
Ih Sa ie al (5.10) 


where k is a suitable scalling factor which can be evaluated on the basis 
of some other measurement (temperature calculated from plasma bulk 
electrical conductivity). Finally, assuming that Ne is already measured 
by some technique, by interferometry say, then 
2 

) 


x 2 
Te Sak No /Un 


: (5.11) 


This latter formulae will be used in Section 6.2.4. 


5.3.9 Spot-Track-Pattern 


The unipolar pulsed discharges leave spot track patterns on the 
surface of cathode and anode. Such pattern can be used as a source of 
information regarding discharge dynamics (Kvartskhava et al. 1966). 

Investigation of these effects has the following steps: 

i) the accelerator is disassembled and the electrodes are polished, 

ii) the accelerator is assembled and fired n-times under the same 
conditions, and 

iii) then the gun is disassembled and enlarged pictures of the spot- 


tracks are taken. 
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A Graflec Camera and type 57 Polaroid film (5 x 4) are used for 
investigation of this effect as elaborated in Section 6.2.1 (Figs. 


Ot 47 and *6.075) 


Note that this method is very simple and cheap although rarely used. 
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CHAPTER VI RESULTS OF THE MEASUREMENTS 


A shock-tube model was formulated in Chapter 3, and its quasi-one- 
dimensional version solved numerically in Chapter 4. Reasonable agreement 
has been found between those calculations and detailed measurements of 
Lie et al. (1967). The aim of the experimental part was to perform 
measurements similar to Lie but make the comparison in more detail. 
During the course of the experimental work it has been found that the 
behaviour of our device is dominated by completely different mechanisms. 
Because investigation of these phenomena are of increased interest in 
current world-wide experimental effort, we decided to pursue these new 
aspects. The importance of this (insulator ablation) is given by the 
application of such devices as satellite thrustors and by the absence of 
a valid theoretical model. 

Section 6.1 describes some details of the discharge-initiation 
regime in our device. Sections 6.2 and 6.3 are devoted to experiments 


with the CPA having two insulator geometries. 
6.1 DISCHARGE-INITIATION REGIME 


A general description of the discharge-initiation regime is given in 
Sections 2.1.1 and 2.2.1. For our arrangement, the "plasmoid" consists 
of the ionized ablated teflon, ionized helium, and ionized impurities 


from walls, gas, electrodes etc. 
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Schematic diagram of the apparatus for initial breakdown. 
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Fig.6.2 Verification of the preionization circuit;oscillograms of the HV-probe signals 


from four consecutive runs.HV:2 kV/div,magnetic probe:2 V/div,l usec/div,axial location: 
2 cm. 
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A simplified schematic of the apparatus, valid for the first 
microsecond is given in Fig. 6.1. One should mention the importance of 
the 100 ke resistors for this regime because they are connected in 
parallel to the teflon-surface-breakdown path and improve the jitter in 
breakdown time. The breakdown sequence is given by #1, #2, #3 as shown 
in Fig. 6.1. The original purpose of the trigger pins in the outer 
electrode was the creation of a current layer at this POSITION aaalNis, 
“horace. wdid.nat happen. Note that this unsuccessful effort is quite 
common for this particular problem (Popovic, 1969). On the other hand, 
Fig. 6.2 shows that the shot to shot reproducibility has been improved, 
even in our "poor" vacuum conditions. Many authors claim that high-_ 


purity gas conditions are required for good reproducibility (Pert, 1970b). 


6.2 ACCELERATOR WITH INSULATOR I 


Workman (1965) has treated analytically the insulator ablation model 

proposed by Keck (1964) and has shown that the ablation loss can be 
expressed with a single correlation parameter (see also Section 2.4). 
Lie (1970) also found a standing arc near the insulator surface (basic 
assumption of Workman's model), but after investigation the behaviour of 
the CPA having Vycor, Pyrex, Lucite, and Mylar insulators he concluded 
that further study was needed to develop a quantitative theory which 
combines the appropriate snow-plow model with the ablation model of 
Workman. 

In the series of insulator types (glass+plastics), teflon is the most 


ablative material, but has not been systematically investigated. Such 
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detailed measurements are described in this Section, 
6.2.1 Current Distribution 


The slug and snow-plow models of the current sheet motion assume 
that the sheet is very narrow. Figure 6.3 shows the current distribution 
inside the CPA having Insulator I. Evidently completely different 
current distribution is obtained. The standing arc yields insulator 
ablation during the whole discharge time. Typical oscillograms of the 
differential magnetic probe and total current at an axial position of 
22 cm are shown in Fig. 6.3. 

From this oscillogram and others, taken at different axial positions 
we can calculate the current distribution along the electrodes for several 
different times. Results of such calculations are sketched in Fig. 6.4, 
where dI/dx is the current per unit length. The total current I is then 
given by 


60 cm 
ia ii (d1/dx)dx & 70 KA (6.1) 
0 


The latter relation has been used for calibration of the differential 
magnetic probe. 

Note that such a current distribution is completely dominated by the 
insulator ablation and only the term current front has any physical 
meaning. This type of behaviour has been reported by Belan et al. (1970) 


and Kalygin et aZ. (1971) but is "unknown" to most recent papers dealing 
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ee Current distribution for CPA with insulator I 35 Usec/div,V_ = -10 kV,axial position: 
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Fig.6.4 Typical time development of the discharge current inside CPA with Insulator I; 
Ppa) LORE. 
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Fig.6.5 Current front position for CPA with Insulator I;p, = 0.02 to 2.0 Torr. 
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with the insulator ablation effect (Lie, 1970); Vondra et al. 1970). 
None of the latter has measured this "local current distribution" (Fig. 


Gea). 
6.2.2. Current-front Position 


The velocity of the front of the distributed current can be easily 
measured by nonintegrated magnetic probes (Fig. 5.12), or assuming that 
this current front is identical with the luminosity front, one can use a 
photodiode (Fig. 5.10). The latter method has been used for the CPA 
having Insulator I and the result is shown as a function of time and 
pressure in Fig. 6.5. An important conclusion of this measurement is 
that a constant velocity of % 4 cm/usec does not depend on the helium 
pressure in the region .02 to 2 Torr and this can be explained simply 
by assuming that the insulator ablation process supplies a sufficient 


amount of working fluid. 


6.2.3 Electron Density and Temperature 


Electron density and temperature have been measured for CPA having 
Insulator I only using i r diagnostics. (Measurement of the case with 
Insulator II has an increased requirement for noise reduction and speed 
of detection). Figure 6.6 shows a typical oscillogram of the local 
current (diff-magnetic probe) and the fringe beats, measured with the 
Mach-Zehnder interferometer (Fig. 5.14) for 4 = 10.6 um and a plasma 
thickness Ag = 4.2 cm. 

The local currents can be calibrated by using Eq. (6.1). The 


electron density, calculated according to Fig. 5.15, is also shown in 
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Fig.6.6 Typical traces of the local current and fringe beats and electron density (evaluated 
according Fig.5.15)3p, = 1 Torr,V, = -10 kV,axial position:30 cm. 
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Fig.6.7 Current and voltage measurements;5 usec/div,V, = -10 kV. 


Fig. 6.6. This electron density has good time correlation with the local 
current measured at the same axial position. 
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The maximum average electron density is % 1.7 x 10/6 cm 
insufficient to produce measurable absorption (see Section 5.3.4). This 
conclusion was verified experimentally. The average (bulk) electron 


temperature can be evaluated from the average (bulk) electrical conductivity 
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o. For this purpose current and voltage oscillograms have been taken, and 
are shown in Fig. 6.7. From these oscillograms and Fig. 6.5, where the 
Current front position & is shown we can calculate the bulk electrical 
conductivity in the middle layer of the CPA discharge. To approximately 


Calculate the average conductivity for t = 10 usec, read values 


Me ecU0V weel ev ULRA.) 2 4 G0ecm 
from which the calculated o is % 1.5 qv! cm! , Using Spitzer's relation 


for conductivity (Spitzer, 1968): 


ingore 
DOG) Resear ae (6.2) 
Beek ex LT 


where In A 7.0, o is in qr! emz! , and electron temperature iP iSatie. hk. 


From Eq. (6.2) electron temperature at 10 usec, in the middle layer of 


CPA is 
by pecan torx 10°°K (6.3) 


which is an average value for a 40 cm long distributed current plasma. 


The ir bremsstrahlung absorption is negligible but the emission 
is easily detectable. Figure 6.8 shows an oscillogram of the detected 
radiation and calculated electron temperature using Eq. (5.11) and where 
the scaling factor has been determined such that the averaged ifs should 
be equal to £Tg> determined from the conductivity measurement. The 


resulting maximum temperature is % 0.6 eV which is reasonable. 
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Fig.6.8 Oscillograms Of local current,detected ir continuum radiation,and calculated electron 
temperature for known electron density (Fig.5.15) and known average temperature ot 


6.3 . ACCELERATOR WITH INSULATOR II 


The original aim of the experiments was to get a shock-tube behaviour, 
i.e., Create a narrow current sheet. This can be approached by suppress- 
ing the standing arc at the insulator surface. Insulator II of Table 3 
has such a suitable shape. Section 6.3 is devoted to the results of the 
experimental investigation of CPA having Insulator II. 

One "old" but very effective method of investigation has been 
employed: the spot-track-pattern method, described in detail in Section 
5.3.5 which is treated in Section 6.3.1. Only this simple static method 
has shown the geometry of the discharge-initiation regime, acceleration 
regime, and die-out (of the current sheet) regime. 

Current distribution and current-front position are described in 


Sections 6.3.2 and 6.3.3. respectively. 
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6.3.1 Electrode Erosion 


The electrode erosion mechanism is very important for the end-type 


PPA (Section 2.3.3), since the eroded material is the working medium. 

For coaxial and parallel-plate accelerators, however, this phenomenon 

is not desirable. Generally, for plasma accelerators having long electrodes 
and small inter-electrode gap, the plasma-metal interaction plays an 
important role. Tilt of the current sheet in a parallel-plate accelerator 
could possibly be explained in this way (Kvartskhava e¢ al. 1967 ).. s.For 

a CPA Having larger inter-electrode gap (our case), this phenomenon is 


over-ridden by the radial dependence ot the magnetic pressure. 


CURRENT f INSULATOR CENTRAL ELECTRODE 
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Fig.6.9 Region definition sketch for static photography of the spot-tracks on the central 
electrode (Fig.6.10) for p. = 0.6 Torr;oscillogram of the discharge current for ies = 10 kV. 
Region A corresponds to inftial breakdown near the teflon insulator,region B to current 
sheet acceleration,and region C to stopping of the current sheet. 
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The inoving discharge leaves certain characteristic spot-track- 
patterns. Assuming a unipolar current pulse generator is available, and 
"many" discharges are repeated, then the method described in Section 5.3.5 
may be used. 

This type of investigation has been performed for the CPA using 
Insulator II (see Table 3) and a static helium pressure of 0.6 T. 

The anode-spot-track pattern developed on the positive electrode 
after 20 discharges, is shown in Fig. 6.10. This pattern features three 
typical regions along the electrode. The regions denoted A,B,C in Fig. 
6.9, can be identified as: 

(A) initiation region, 

(B) acceleration region, 

(C) die-out region. 

Note that Fig. 6.10 shows only a small part of each region. The Shape of 
Insulator II is such that the standing arc hear its surface is suppressed, 
i.e., region A shows the steady discharge before the acceleration could 
develop. Note that the trigger pins (Fig. 6.1), placed 6 cm from the 


insulator, have no influence on the region A position. The pattern of the 


middle region (B) shows that there is negligible anode erosion during 
acceleration. The pattern of the last region (C) proves that, under 
given conditions, the externally sustained plasma motion stops inside the 
gun (region C has deep anode spots of steady arc structure). The inertial- 
plasmoid motion could be expected to drive a shock wave. 

The cathode-spot-track pattern, developed on the negative electrode 
after 100 discharges, is shown in Fig. 6.11, again for regions A,B, and 
C. The difference is in the shape of the individual anodic and cathodic 


spots. Note that acceleration region B has deep ~ 1 cm long cathode- 
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tracks but almost invisible anode-tracks. The discontinuous form of 
these patterns has been known for a long time (Kvartskhava et al. 1966) 
and is in discrepancy with the usual "smooth" current sheet assumed in 


most available models. 
Grose “Current: Distribution 


The technique used is the same as in Section 6.2.1. The CPA having 
Insulator II has narrower current sheet than that having Insulator I, but 
still it is not the required few millimeters thickness of the snow-plow 
models. Figure 6.12 shows measured local current distribution measured 


at a fixed axial position 2 = 10 cm for helium pressures of 1, 2, and 


4 Torr. Because the measured velocity is = 1 cm/usec for a pressure of 


1 Torr (Section 6.3.3) then the layer thickness is roughly 2 cm. 
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100 kA/div 2 
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Fig.6.12 Current sheet distribution for CPA with Insulator II; 2 usec/em,V, = +10 kV. 
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6.3.3. Current front Position 


The technique used for the front position determination is that 
illustrated in Fig. 5.10, i.e., the method uses one photodiode, two 
irises, and multiexposure of oscillograms. Results are in fl GO seals 
In comparison to the CPA having Insulator I, these results show that gas 
pressure has a weak influence on the current front position, but the 


| snow-plow velocities (% 6 cm/usec) are still not reached. 
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Fig.6.13 Current front position of CPA with Insulator II;V, = +10 kV. 


fut 





soit et no teanliemeed nehtireg Joo att 30% Beapeaupiad ised WAY” Pre 
ow? ,aboTberudg sna seen Godddn og ,.a/? (OR «prt at badeaseuf ty 
f,8 QT ot ove edfucell .ateigol i rose To s7emues Ton Ris geeePr 
hy Set) Wide 2iueed ceete ,1 Jose tude! tort by tel NSS ort? a nie t/seyinory nt 
of thd .nots eq snow? tienda! 8 eae anew 6 2a sweesng oly 


% 


HA IES ror ' in Gh { cbt Vere a #) 25) itaolsy wolg-wore 





: eiete ia, | : 
’ Laie? 74 7 ' 
: 
| rae . t a 
‘ 4 “ Me ~*~ \ : i 
f ~~ Mas x ‘ a 
; r ™ . ~~ oo ™/~ * -s 7 } 
~, A ’ ~ : 
. ' A. , ~ q i a > 
\ ee 
$ 
; on 
i - 
Pe ee eects — a ll - 
te te gi re 
ee a eT | a pat ee 
: ; . ic a 
os i? ‘ 7 
‘ 
| acne is enema, ‘a aya 
7 - ~ : » a ny 7 
La i ips rs en : rs rc 





oy asta a wr Deity or 
7 : .) ; ah . 


if 
yy 8 
ee 
; 


102 


CHAPTER VII RESUME AND CONCLUSIONS 


Coaxial plasma accelerators have been described by introducing a 
set of regimes (kind of simple behaviour), regime models, and types 
(kind of accelerators). There is no such complete review available in 
the present literature. 

A more detailed theoretical description was given to one type of 
formulation of the axisymmetric-steady-shock-tube regime, and two 
versions [Y(r,z) and Y(r,2(r)] were developed. 

The two-dimensional model, Y(r,z), was approached in a very 
simplified version because of the complexity of the numerical methods 
involved. Verification of the resulting procedure was obtained but was 
found to be rather sophisticated, expensive, and parameter sensitive for 
the original purpose. 

The quasi-one-dimensional model, Y(r,z(r)), leads to a simultaneous 
set of ordinary differential equations and nonlinear equations. Numerical 
procedures involved were found to be relatively simple and the solutions 
for helium (1 Torr, velocities 1 to 10 cm/ysec) was given. Comparison with 
available experimental data is reasonable. 

Accelerators with two types of teflon insulator were investigated 
experimentally. The first type features a standing arc near its surface 
and the resulting acceleration regime is completely dominated by the 
insulator ablation process. Related velocities, electron densities and 


temperatures, and currents are measured. The second type of insulator 
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has a special shape for standing arc suppression and the related 
acceleration regime was found to be midway between the insulator-ablation 
and shock-tube regimes. The pure shock-tube regime has not been reached. 
On the basis of this and other recent experimental work, the only 
available insulator ablation model (Workman,1965) should be extended. 
A unipolar current pulse generator (100 kA, 35 usec), developed 


for this study is of general applicability in pulsed plasma investigations. 
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